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ABSTRACT 


The  feasibility  of  stabilizing  gas  bearing-rotor  system  by  electromagnetic 
me  an  o  i.ovc  a  t  lcu  ana  iy  t  lc  a  1  ly .  Two  devices  appeared  feasible,  namely, 

a  unidirectional  device  producing  a  magnetic  force  to  load  the  bearing  with¬ 
out  increasing  the  rotor  mass,  and  an  active  device  producing  a  controlled 
electromagnetic  force  always  opposing  the  motion  of  the  shaft.  A  numerical 
example  shows  that  the  power  loss  from  using  either  device  together' with  plain 
journal  bearings  compares  favorably  with  tilting-pad  gas  bearings- 


(This  abstract  is  subject  to  special  export  controls 
and  each  transmittal  to  foreign  governments  or 
foreign  nationals  may  be  made  only  with  prior  approval 
of  the  Air  Force  Aero  Propulsion  Laboratory  (APFL) , 
Wright-Patterson  Air  Force  Base,  Ohio  45433.) 
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SECTION  I 


INTRODUCTION 


In  high  speed  machinery  applications,  gas  bearings  have  many  advantages  such 
as  low  friction,  high  stiffness,  insensitivity  to  extreme  thermal  environment, 
radioactive  ’'contamination",  and  so  on.  However,  a  gas  bearing,  because  of 
its  low  damping  characteristic,  is  susceptive  to  self-excited  vibration  which 
appears  in  the  form  of  journal  whirl  motion.  This  motion,  once  started,  will 
grow  in  amplitude  until  it  either  reaches  a  finite  amplitude  or  results  in 
solid  contact  between  the  journal  and  the  bearing  and  thus  causing  bearing 
failure.  Even  the  finite  amplitude  whirl  may  sometimes  be  considered  unsat¬ 
isfactory  depending  cn  the  particular  application  and  the  magnitude  of  the 
whirl  amplitude. 

Shallow  spiral  grooves  engraved  on  either  the  journal  or  the  bearing  surface 
were  found  both  theoretically  and  experimentally  to  be  effective  in  reducing 
whirl  instability  [Ref.  1,  2  and  3l*.  However,  in  some  range  of  operation, 
even  a  grooved  bearing  can  have  whirl  instability.  Tilting-pad  gas  bearings 
have  been  used  successfully  in  many  applications.  Although  tilting-pad  gas 
bearings  can  be  designed  to  operate  quite  stably,  the  designs  are  complex  and 
power  loss  is  high.  Therefore,  other  means  of  improving  the  stability  of 
journal  bearings  are  desirable.  One  possibility  is  to  utilize  magnetic  forces. 

In  the  literature,  magnetic  suspension  systems  have  been  investigated  [/*]  for 
apace  applications.  The  advantage  of  a  magnetic  bearing  over  a  gas  bearing  is 
its  ability  to  support  load  even  in  the  absence  of  any  lubricant.  This  makes 
magnetic  bearings  particularly  suitable  for  suspension  systems  in  high  vacuum 
environment.  However,  magnetic  bearings  generally  have  much  lower  load  capac¬ 
ity  and  stiffness  than  gas  bearings. 

In  view  of  the  above  it  appears  logical  to  use  the  advantages  of  both  a  gas 
bearing  and  magnetic  forces,  i.e.,  we  will  use  gas  bearings  to  provide  load 
capacity  and  stiffness,  and"  magnetic  forces,  hopefully  to  supply  the  damping 
for  bearing  stabilization. 

*Numbers  in  brackets  refer  to  references  at  the  end  of  this  report. 
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The  purpose  of  this  effort,  summarized  by  this  report,  was  to  investigate  the 
feasibility  of  applying  electromagnetic  forces  to  stabilize  gas  journal  bear¬ 
ings.  Plain  cylindrical  journal  bearings  were  chosen  ft'  investigation  be¬ 
cause  of  their  simplicity  in  analysis  and  because  thry  served  the  purpose  of 
demonstrating  how  the  whirl  instability  can  be  suppressed  by  electormagnetlc 
devices. 


SSi  MOV  Li 

■THEORETICAL  gACKLROUND 

As  stated  in  the  Introduction-  a  plain  ylcndncal  journal  bearing  was  chosen 
for  investigation.  We  shall  briefly  recaci tuia'e  in  this  section  the  theoret¬ 
ical  background  of  stability  analysis. 

Following  the  approach  [i],  assume  tnut  the  journal  has  a  steady  sf3te  displace 

ment  e  from  concentric  position.  Let  us  choose  the  x-axis  to  coincide  with 
xo 

the  line  connecting  the  concentric  position  and  the  displaced  equilibrium  posi¬ 
tion,  and  the  y-axia,  perpendicular  to  it  (See  Fig.  1).  Let  the  journal  have 
small  oscillations  with  frequency  (1  and  amplitudes  e  and  e  ¥  In  the  x  and  y 
directions  respectively.  Then,  from  Fig.  1..  it  is  illustrated  that  due  to  the 
steady-state  displacement  e  ,  there  is  a  gap  variation  e^o  cos  9  (detailed 
derivation  is  provided  on  Pages  41  and  42  of  Ref.  [10]> .  Similarly,  for  the  x 
and  y  oscillations,  e  .  cos  At  and  e  _  coa  At,  the  gap  variations  are  e  cos  At 
cos  0  and  e  cos  At  sin  0  respectively.  Thus,  the  dimensionless  bearing  gap 
can  be  expressed  by 
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where 


Zxx  -  Uxx  +  1  vxx*  etc- 


(4) 


(5) 


The  Z  matrix  is  the  mechanical  impedance  of  the  bearing;  the  U  and  V 
are  the  direct  dynamic  stiffness  and  damping  in  the  x  direction,  and  z 
and  Zyx  are  cross-coupling  terms.  The  Z:s,  F  and  G  are  given  in 

[5]  and  listed  in  Appendix  I  for  easy  reference;  they  are  functions  of 
€xq,  f  and  a). 


In  general,  the  motion  of  the  Journal  can  be  expressed  by  e  e3T  and  e 
e  ,  with  t  «  cot.  At  the  onset  of  instability,  s  must  be  purely  imaginary, 
because  a  real  part  would  make  an  oscillatory  motion  either  grow  or  decay 
depending  on  the  sign  of  the  real  part.  Thus,  we  set 
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where  is  the  critical  speed  ratio. 
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And  the  equations  of  motion  are 
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Sl.c.  (9)  i.  homogeneous ,  the  dee.rmln.ot  of  the  coefficient  metrl, 
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The  u' b  and  V'i  are  function*  of  f  aa  indicated  in  Appendix  I  and  Ref  [5]. 
For  a  given  rotational  speed  cu  and  steadys^ate  displacement  c  ,,  Eq.  (11) 
determines  the  critical  speed  ratio  f  and  Eq.  (12)  determines  the 
critical  mass  parameter  m  . 


It  is  customary  in  stability  analysis  to  define  a  threshold  speed  from  the 
critical  mass. 


Threshold  speed 
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The  critical  mass  parameter  and  the  threshold  speed  are  the  upper  limits  for 
stable  operation  of  the  bearing.  The  actual  values  should  be  desigr.ti  to 
stay  below  their  upper  limits  to  have  a  stability  margin. 
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In- t.he  "last  section,  we  have  describee!  that  the  stability  of  a  journal 
bearing  can  be  expressed  by  a  critical  epees  and  a  critical  mass.  Ihe 
mass  of  the  rotor  must  be  kept  below  the  critical  mass  for  a  sc-:,  il  a 

operation.  -v 


It  is  shown  in  Ref.  2  that  an  unloaded  plain  cylindrical  journal  bearing 
is  always  unstable  because  it  he-  a  zero*  critical  mass.  This  situation 
is  inevitable  if  a  bearing  is  co  be  operated  in  a  zero-g  environment. 

One  possible  means  to  make  the  bearing  stable  is  to  apply  a  unidirectional 
magnetic  force  so  that  the  bearing  will  be  run  at  an  eccentric  rather  char: 
the  concentric  position.  Detail  consideration  to  this  approath  will  be 
presented  later  in  this  section. 


Although  we  can  calculate  for  a  leaded  journal  bearing  (grooved  or  uogrooved) 
the  critical  mass  for  instability,  in  many  cases,  the  critical  mass  turn  a" 
out  to  be  lower  thary  the  rotor  ms  as.  and  the-  system  would  be  unseat  ie . 
Therefore,  it  Is  dedirsble  to  devisa  means  to  increase  the  critical  mas; 
and  thereby  make  the  .system  more  stable.  One  such  device  will  be  shown  to 
be  very  effective  in  achieving  the  above  purpose;  it  will  ba  called  ar. 
active  electromagnetic  device  (see  Fig.  8). 

Electromagnetic  Relatior snips 

Let  us  first  recapi tulate  some  caeic  elertro'magr.etic  formulae.  ~~  at. 
electromagnet,  the  magneto-motive  force  (m.m.f.)  is  equal  to  the  ampere- 
turn  of  the  winding, 
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=  air  gap  area  perpendicular  to  magnetic  flux,  in^ 
l ■  ™  length  of  air  gap,  i-_. 

I k  «  mean  length  of  winding  per  turn,  in. 

A^  *  area  of  conductor,  circultt  mils 
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saturation  flux  density  for  this  ca.it  alloy  is  50?000  lmea/in  .  As  will  bt 

2 

seen  later,  the  actual  flux  density  u  designed  below  40,000  lines/in  . 


Now,  for  the  purpose  of  an  order-of  magnitude  analysis,  suppose  that  a  high 
permeability  material  such  as  supermalloy  is  chosen  for  both  the  core  and  the 
rotor,  we  may  assume 

i)  uc  »  *  140,000  line/in2y^rap. turn/in 

and  ii)  that  the  "flux  passage  length"  through  both  the  core  and  the  totor 
is 

A  A 

*c  A*  +  A  *“  -  20  l“- 

c  r 


The  permeability  of  air  ia  UQ  *  1  Gauas/Oersted  -  2.54  line/amp.turn  in. 
Assuming  the  air  gap  to  be  0.01  in.,  then, 

Ac  20  x  2.54  1  , 

21  *  140,000  x  2  x  0.01  **  56 *  < 

© 


Therefore,  for  this  analysis,  the  combined  reluctance  of  the  core  and  the  rotor 
can  be  neglected  when  compared  to  the  reluctance  of  the  two  sir  gaps.  The  re¬ 
luctance  of  the  magnetic  circuit  is,  therefore, 


Using  the  above  expression  for  the  reluctance  of  a  magnetic  circuit  and  follow¬ 
ing  the  derivation  shown  on  page  33  of  Ref.  9,  we  obtain  (instead  of  Eq.  (14) 
on  page  33  of  Ref.  9)  a  relationship  between  flux  density  and  m.m. f , , 


Ni  -  0,313  B  (2lg)  «  0.626  B  (16) 

Note  that  given  by  Eq.  (15)  ia  the  magnetic  force  produced  by  one  (magnetic 

air)  gap.  The  total  magnetic  force  will  obviously  be  2?  . 

in 
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Since  <:hers  are 


two  journal  bearings  -  one  at  each  enc  of  the  shaft  -  each  beating  will  be 

forced  to  carry  a  magnetic  snatc  loaaing  o£  5  1-.  audit! or.  tc  other 

m 

shaft  loads. 

1.  A  Device  Producing  a  Unidirectional  Magnetic  Force 

A  schematic  diagram  of  this  device  is  shown  In  Fig.  2.  The  magnetic  force  is 
produced  by  a  d-c  electromagnet.  Magnetic  lines  are  emitted  from  one  "leg" 
of  the  magnet  (north  pole)  to  the  rotor  and  bac a  into  the  other  "leg"  ( south 
pole).  The  rotor  and  the  magnet  should-,  of  course,  be  made  of  magnetically 

permeable  material. 


Under  steady-state  condition,  the  current  in  the  winding  is 
E 

‘o  •  li2 


(17) 


where  Eq  ■  steady-state  electric  potential  [volt] 
R  *  resistance,  [ohm] 


Knowing  1,  one  can  easily  calculate  B  from  Eq.  (16)  and  F  from  Eq.  (15).  Note 

m 

that  the  magnetic  force  F^  is  always  an  attractive  force.  A  magnetic  stiffness 
can  be  defined  as  the  increment  of  magnetic  repulsive  force  per  unit,  decrement 
of  the  air  gap. 

dF 


ra 


d(-F  ) 
in 

-dl 


m 


d/, 


( 18) 


8 


Using  Eqs.  (15)  and  (16)  and  carrying  out  the  differentiation,  it  is  readily 
obtained  that 
4F 


m 


m 


2 1 


( 19'; 


8 


Note,  first  of  all,  that  the  magnetic  stiffness  is  negative.  This  is  an  un¬ 
desirable  feature  because  the  stiffness  of  the  rotor-bearing  system  will  be 
decreased  by  the  amount  ^'tr/2-tg* 
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In  order  to  keep  the  (regritudr.  of  4^  sirs'.  1,  it  is  desiraoie  to  have  3  -c  -  -  i  1  <- 

targe  i  ;  this,  of  course,  requires  a  proportionally  iarg;  Nt  to  driv/C-  the  same 

6 

amount  of  magnetic  flux  through  the  a<r  gsp  It  ir  to  be  noted  that  k  is  the 
static  magnetic  stiffness.  Lndar  dynamic  condition  vi.e.,  the  gap  has  a  sinu¬ 
soidal  variation  with  time},  the  situation  is  different.  First  of  all,  the 
electric  current  is  no  longer  E/R  1.-st  us  choose  a  coordinate  syscem  £  as  shown 
in  Fig.  1.  Assume  that  the  motion  of  the  rotor  ir.  the  £-direttior  be  represented 


by  e cos  fit  with  e_  much  smaller  than  £  . 

ft  g 


i  +  l'  \  V  «  i 

go  g  g  go 


Then  the  air  gap  can  be  expressed  by 

(20) 


where 


*g  "  "  e5t  cos  nt 


(21) 


The  inductance  of  the  electric  circuit  is  given  by, 
L  «  N  “  10  ^  henry 


(22  > 


But, 


CP 


Thus, 


BA  » 


Ni 


g  0, 313(2 £g)  g 


•  3i 


3.19  x  10 


x  n2a 

-8 _ & 


21 


(23) 

Note  that  L  is  time  dependent.  Therefore,  instead  of  the  usual  current  equation 


g 


L  ~  +  R  i  *  E 
at  o 


we  should  use 


ft  <Li)  +  R  1  -  Eo 


(24) 
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Using  Equations  (20'  and  (23 ) Eq.  '■ ?*•>  c  en  Sc  lieu  by  neglecting  quad¬ 

ratic  and  higher  order  effects  of  the  air  gap  perturbation.  First,  Eq.  (23) 
can  be  written  as 


L  «  L  +  L' 
o 


(25) 


where 


n2a 


L  »  3.19  x  10 
o 


L'  -  L 

lgo  0 


21 


g° 


(26) 


Since  the  current  is  alao  varying  with  time,  we  can  write 


i  -  i  +  i  V 
o 


(27) 


where  i  ia  the  steady-state  current  and  1',  the  perturbation.  Now  Eq.  (24) 
o 

becomes 


fr  (L  +  L')  (i  +  i')  +  R  (i  +  i’) 
at  o  o  O 


(28) 


which  can  be  separated  into  a  steady- state  equation, 


R  i  «  E 
o  o 


(29) 


and  a  perturbation  equation, 


1  4—"  +  i  4~  +  R  i 1  »  0 

o  at  o  ai 


^  30) 


dl *  dL* 

Note  that  in  Eq.  (30),  we  have  neglected  L1  and  i1  ,  because  these 
two  terms  are  small  when  compared  respectively  with  the  first  two  terms  in 
I,.  (30),  i.e.  snd  1 1  <  <  i„  f1- 
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f 


Equation  (29)  indicates  that  the  steady- s: ate  current  is  ~  which  agrees  with 

K 

Eq.  (17)  and  is  obviously  valid.  Equation  (30)  is  the  equation  for  the 
pcrturb'atiun  cut  rent.  Rearranging  (30)  .  we  obtain 


di  ,  R 


—  i'  *  -  3~ 

L  "o  dt  ' 


Corresponding  to  a  periodic  perturbed  motion  of  the  rotor  sucn  that 


l  “  l  -e  cos  n  t 
g  go 


the  solution  of  Eq.  (31)  is  [Ref.  12] 


i’  "  *  7^  0-  "2 - 2 - 2  "in  0  C 

°  4go  |_n  >^R/Lo)2 


n2  +  (r/l  )2 

o 


cos  0  t 


+  i"  exp  (-  j-  t)  i  32 

The  transient  term,  which  decays  exponentially,  depends  on  the  precise  initial 

egt 

condition;  for  Instance,  assuming  i'  (t«0)  -  0,  we  have  i"  ■  l  — * —  x 

n  °  *8° 

-  -**-  .  — ,  t  in  any  case,  the  time  constant,  L  /R,  in  a  typical  design. 

a  +  (r/lo)2  0 

would  be  very  small  as  will  be  illustrated  in  a  numerical  example  on  pages  1' 
and  13. 


From  Eq.  ( 16) 


Ni  •  ^o*1') 

0.626  i  “  0.626  (/  +/') 

*g  v*go  *g' 


»o  <  1  *f 

o 


where 


■Ni 


c  0.626  i 


34, 


'go 


Finally,  we  can  compute  F  from  i,q.  (13‘,  and  by  neglecting  terms  involving  pro- 

m 

ducts  of  perturbation  quantities. 


Ag. 

.  ,6 

72s  i0 


i  *b 

'  1  +  2  --  ■  2X1 


*go 


Let  F  *  F  +  F  1 
oi  mo  m 


(35) 


Then, 


mo 


Bo  A* 

o  g 

72x10* 


(35*) 


F  1 


F™,(2f 

o 


lo 

*  2 

go 


(35b) 


Substituting  Eqs.  (21)  and  (32;  into  Eq.  ( 3 5 b >  and  assuming  that  the  transient 
term  has  already  become  negligible.,  we  iotain 


F‘  -  F  2  ~  «... 
®  mo  lgQ  It 


fJR/L 


2  2 
n  +  (R/Lp) 


-(R/L 


2 


n2  +  (r/lo)2 


sin  CX 


cos  n  t 


(36) 


where  « 


C 


.  ^ 

c 


dynamic  eccentricity  ratio 


radial  bearing  clearance 


F^  is  the  perturbation  force  due  tc  a  shaft  oscillation, -e gt  cos  fit.  Thus, 

the  in-phase  and  out-of-phase  components  of  F'  represent  respectively  the 

m 

dynamic  stiffness  and  damping  induced  magnetically. 
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In  order  to  express  the  magnetic  force  In  terms  of  stiffness  and  dampine.  we 
write  Eq.  (4)  in  che  §T)  coordinates. 


J_  pv* 

11 ».  [<V._ 


Z*!5  Z"5T1  ["‘St 

Z!HT15  ZfnlwJ  I  *T|t 


where  Z  «  U  +  i  V 

m'5  mS?  mC5  >  etc*  The  subscript,  m9  indicat. 
due  to  magnetic  forces.  From  Eq.  (37)  we  have 


es  quantities 


ttDL  p 


•[v  '"»«  a  1  -  'V«  «”  n  ' 


+  >  <U.,T,  "■  n  '  -  vm{1 .1" 


n  t  )l 


t£  ..  Identify  of  Eq.  (36)  „lth  (I  ^  of  Eq.  (38).  .„d  equ.te.te™,  of  like 


combinations  of  *  cos  fit  etc.,  we  obtain 


U  «  -  2  F  ——  ^ 

m5S  mo  ig0  x  + 


c  (flL  /R)  • 

V  -  2  F  - 2 - 

"  480  1  +  (flLo/R)  2 


where 


roo  ttDL  p 


W2  +  W2 
XO  yo 


? 


'  * *•*  4*«ofc. «« . 


Iii  nicer  co  express  the  magnetic  force  in  terms  of  stiffness  and  damping,  we 


Observe  chat  as  0  -»  0  -'of  Q<  <  tnr  magnetic  dynamic  stiffness 

approaches  the  steady -state  stif.fr, is;,  a:-  .lies-  by  Eq.  il3)  except  a  constant 
factor  due  to  nou-dimensior.siizscior .  ard  :he  magnetic  dynamic  damning 
approaches  to  zero.  Therefore,  cne  qja:-i  static  resist;,  as  expected, 
agree  with  the  static  results  ob‘aired  octore. 

Eolation  (39)  indicates  chat  t.he  electromagnet  aside  from  providing  a 
lo  Hng  on  the  bearing  and  a  negative  stiffness,  yields  a  positive  damping 
dynamically;  both  are  inversely  proportional  to  the  magnetic  gap.  They  are 

a* 

plotted  against  fl  L  /R  in  fig-  3.  Since  the  damping  helps  to  stabilize 
whereas  the  negative  stiffness  tends  to  deste.oi.ize,  it  is  desirable  to  have 

large  values  of  V  __  but  keep  the  magnitude  of  __  small.  i.n  order  to  do 

®55  ®55 

so,  the  parameter  0  L  /R  should  be  designed  to  be  greater  than  unity  as 

o 

can  be  seen  from  Fig  3 . 

Recall  that  g  Is  the  direction  of  magnetic  loading.  If  the  m-rgnet  fare  is 
designed  to  have  a  small  wrap  angle,  so  that  Che  reluctance  is  not 
appreciably  changed  with  s  shaft  motion  ir  the  Tl~ditection  .  vs  have 


'Jm7]T]  "  Vum  “  C  >C) 

All  the  cross-coupling  terms  are  also  zero. 

Now  we  can  express  the  magnetic  dynamic  stiffness  and  damping  in  the  xy 
coordinates. 


{  Umxx,  Vmxx  }  “  {  ^mgg,  Vmgg  }  cos2“ 
{  Umyy,  Vmyy  }  *  {°mgg;  ’'mgg  }  sln  a 


{  Umxy,  Vmxy}  *  VJ'mg§,  ''mgg  } 
l^myx ,  ^myxi 


mgg ,  y«K  ; sln  *  cos  * 
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The  angle  3  is  Che  angle  of  rotation  from  the  gt]- coordinates  to  the  xy- 
coordinates  as  shown  in  Fig.  1 


The  dynamic  stiffness  and  damping  of  the  rotor- bearing  system  iB  equal  to  the 
sum  of  the  corresponding  terms  of  the  bearing  (U  ,  •••)  and  the  electro¬ 
magnet  (U  ,  . ..;. 
mxx 


(U  )  -  U  +  U 

xx  i  xx  mxx 


(U  )  «  U  +  U  ,  etc. 

xy  a  xy  mxy’ 


(42) 


where  the  subscript  s  indicates  the  system.  Knowing  the  system  dynamic  stiff¬ 
ness  and  damping,  the  critical  speed  ratio  and  the  critical  mass  can  be  calcu¬ 
lated  from  Eqs.  (11)  and  (12). 


Example  —  Plain  journal  bearing  In  zero-g  field. 


Input  D  -  1  in. 

Pa  *  14.7  psi 

/  -  0.1  in. 

go 

C  »  0.001  in. 
u  *  0.27  x  10"8  lb  sec/in2 
A  -  0.1 


A  *  0.2  in2 

S 

A^  *  1000  cir.  mil. 

jtk  «  2.5  in. 

i  ■  1  amp 
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From  A 


6  u  co  .  D  .  ‘ 
pa  (2C> 


u,  c 


we  calculate  a  *  363  rad/at'., 


To  displace  this  bearing  to  an  eccentricity  of  cXQ  "  0.2,  we  need  a  bearing 


force  of,  from  Eq.  (3), 


W 


xo 


W 


Then 


yo 

a 

F 

mo 


-  0.696  x  10 
0.24  x  10'2 
88.3° 

“v/w2  +  w2 

V  xo  yo 


-4 


-  0.0024 


which  corresponds  to  a  dimensional  loading  of  F  ■  F  ^  ’  ttDL  p  »  0.111  lb. 

■  mo  mo  r  a 


From  Eq.  (15)  and  Eq.  (16)  we  calculate 

B  -  6326  lines/ln2 
N  -  392  turns 

The  inductance  is,  from  Eq.  (26) 


M2A 


“  -  ,  7s—  3.19  x  10*8  *=  0.496  x  10-2  henry 

°  ho 

If  copper  is  used  for  the  coil  winding,  the  resistance  is 
R  -  -  0.85  ohm 

Thus,  the  transient  time  constant  is 


Lo  _  0.496  x  10~2 


R 


0.85 


-  0.00583  sec. 


snd  it  would  take  only  0.027  sec.  for  the  magnitude  of  the  transient  term  in 
Eq.  (32)  to  be  reduced  to  1/100  of  its  initial  value. 
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<:  — - * -**  -  <«> 

actually  ah.  ""*  °y  f  *  »•««  (which  i. 

y  he  critical  aped  ratio),  than,  Li  v  can  be  calculate,! 

'**  E’-  <”>  «*  »,*,  etc.,  fro.  Eq.  (41).  55  "« 


! 

f 

l 

i  u 

i  mxy 

u 

mxx 

*  -0. 19  x  10'7 

U 

myx 

*  —0.656  x  10“  ® 

% 

V 

myy 

-  -0.227  x  10*4 

i 

V 

mxx 

*  0.201  x  10“ 7 

V 

S  ““y 

V 

myx 

■  0.695  x  10“ 6 

1 

1 

V 

myy 

”  0.24  x  1 0”4 

r  2  Ch‘  “  b'  “lC0Ut'‘l  *>•  <*»  ...-in, 

>:x  on  have  already  been  calculated  from  Ref.  [5], 


<Uxx>a 

-  0.716  x  io*3 

(Uxy>a 

-  0.120  x  10“ 1 

<V- 

“  -0.123  x  10’ 1 

<0  ) 
yy  * 

-  0.664  x  10“ 3 

(v  ) 

xx'a 

■  0.122  x  10* 1 

(V  ) 
xy  a 

•  -0.693  x  10“ 3 

(V  ) 
yx  e 

*  0.681  x  10‘3 

(V  ) 

yy  « 

■  0.120  x  10“l 

From  Eq.  (13) 


Z,  - 


-0.586  x  10" 3 
0.667  x  10*4 
— 0. 1374  x  10’ 2 
■“0.2425  x  10* 1 
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Substituting  into  Eq  ,11';  .  we  have  V  *  '  ■*  Z  ~  Q.\?t  x 

°  xx  s  yy  s  1 


Thus,  the  vain*  of  f  *  0.4995  is  inde.«sd  he  critical  speed  racit 
less  critical  mass  is  from  Eq.  ;12) 


r>p  di  ncnsion- 


m  -  0.108  x  10 

c 


Threshold  speed  -  0.067 


For  a  plain  journal  bearing  with  L/D  »  1,  C  ~  0.001  in  and  e^o  »■  0.2,  the  thres¬ 
hold  speed  is  plotted  against  A  for  various  values  of  A  ^  in  Fig.  4.  It  is  seen 
that  the  threshold  speed  increases  with  t  in  the  small  A  region.  When  l  -*  to . 
the  stability  curve  becomes  that  of  a  gravitationally  loaded  bearing,  the  number 
of  ampere-turns  required  to  magnetically  load  the  bearing  is  of  course  infinite 

(because  A  m  w) .  For  moderate  and  high  A  the  magnetically  loaded  bearing  is  as 
go 

good  as  the  gravity  loaded  bearing.  Similar  stability  maps  for  -0.4  and  06 
are  shown  in  Figs.  5  and  6.,  respectively.  In  Fig.  7  the  number  *r  terns,  N  is 
plotted  against  A  for  the  same  Journal  bearing  at  e  ^  *  0  6,  the  electric  current 
in  the  coil  is  assumed  to  be  fixed  at  one  amp. 


2.  An  Active  Electromagnetic  Device 

In  an  active  electromagnetic  device,  the  motion  of  the  shaft  is  sensed  by  capac 
itance  probes.  (See  Fig.  8).  Since  the  whirl  motion  is  two-dimensional ,  two 
probes  are  needed  and  should-be- placed.  90  degrees  from  each  other  The  output 
of  the  sensing  probes  will  be  amplified,  then  sent  from  the  amplifiers  :o  the 
windings  of  the  electromagnets  placed  around  the  shaft.  The  electrical  system 
will  be  connected  in  such  a  way  that  the  electromagnets  always  exert  forces 
opposing  the  motion  of  the  shaft.  Since  an  electromagnet  can  only  exert  attrac¬ 
tive  forces  on  the  shaft,  two  electromagnets  are  needed  in  tne  horizonra-.  dire.c 
tion  and  two  in  the  vertical  direction.  Therefore,  a  total  of  four  ‘lectromag- 
nets  are  required  for  the  device.  Figure  8  illustrates  a  schematic  diagram  of 
the  device.  It  is  seen  that  there  are  two  independent,  but  identical  subsystems 
-  -  one  to  control  the  x-motion,  and  the  other  to  control  the  y-tnotion  Each 
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i 


subsystem  consists  or  or  =  capa-oit  at  *e  prooe..  one  di  frerentiatot  .  ;r.e  ^mpl; - 
fier,  one  signal  apU».n*r  4-4  two  0  lc:rro?.->g..c:s.  It  is  necessary  no  ;.«•/**• 
tigate  the  chararteri scics  of  or.^  -i.*b system  only.  Let  the  di apt a cement  in 
the  x-dlrection  be  represented  by 

e  *  e  +  e  cos  fit  ''63)' 

x  xo  xr. 

where  e  Is  the  stssa;-  state  di spl*-  ament  determined  by  the  loading  and  the 
xo  r  J  -o 

bearing  characteristics,  and  e  is  the  peak  amplitude  of  the  time  dependent 
term.  Because  the  loading  may  change  either  in  magnitude  or  in  direction  or 
both,  exQ  would  change  accordingly.  But  as  far  as  improving  the  system  sta¬ 
bility  is  concerned,  we  have  no  direct  interest  in  e^Q.  As  a  matter  of  fact, 
it  is  desirable  to  sense  only  the  time  dependent  part  of  the  displacement,  to 
that  each  time  there  is  a  new  we  do  not  have  to  readjust  the'system. 

However,  a  displacement  probe  can  only  read  e  .  This  is  why  a  differentiator 

•<  X 

is  needed;  it  differentiates  with  respect  to  time  the  output  of  the 

m 

displacement  probe.  The  input  to  the  amplifier  is,  therefore,  or  -  eyt  fi 
sin  Ot.  Finally,  the  algnal  splitter  is  to  decide  which  electromagnet;  to 
energize  with  the  output  of  the  amplifier.  Physically,  it  is  clear  that  or 
ex  is  negative,  electromagnet  CD  should  be  energized,  or  electromagnet  (3) 
if  otherwise.  A  description  on  the  differentiator  and  the  9igi:al  splitter 
is  given  in  Appendix  II. 

The  output  of  the  differentiator  or  che  inpu*  to  the  amplifier  . a  proportional 
to  the  shaft  center  velocity.  -  e  0  sin  flt..  Let  he  peak  amplitude  of  the 
amplifier  output  ;><- 

If  we  assume  the  time  lag  introduced  by  the.  amplifier  (a*  well  as  the  probe 
and  differentiator)  is  negligibly  small,  then  the  amplifier  output:  pore  it. al 
can  be  represented  by  sin  fit..  In  the  sketer.,  one  full  cycle-  of  E^  sin  fit  1 
shown.  The  positive  portion  (solid  lice'  will  be  channelled  to  energize 
EM  CP  whereas  the  negative  portion  (dotted  line)  Is  to  energize  EM  Q) 
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Only  if  inequality  iuS;  is  satisfied  car  .v~  --,-ctric  c  jrrenc  in  SM  ie  built 
up  s;  iiu  oppitcc  iao  te  time  .ag  ard  daapvd  o  -  quickly  wner>  1:  .  s  nvr  reeded  Ihis 
usually  requires  an  additional  external  electric  resistance  in  t»ri?5  with  the 
winding  Having  made  R  sufficiently  lari?.  :o  satisfy  :i8.  £q  {&?»  can  be  ap¬ 
proximated  by 

A  «s  Ej/R  ;  Y  as  0  (49.) 

In  some  applications,  tne  loss  d:s  to  the.  addition  of  externa:  eiec*nca]  resis¬ 
tance  can  be  quite  large,  it  can  ever,  oy  ■  t *•  d.’miratmg  loss  :n  the  D-’aring  system. 
Therefore,  it  is  interesting  ro  consider  the  alternative  of  using  a  transistor  in 
the  circuit  as  shown  in  the  following  diagram  for  EM  .  Identical  systems 
should  be  applied  to  che  other  three  electromagnets. 


i  t- 

i 


The  driving  voltage  from  signal  sclj.ttrr  i  <  ounen ted  to  pass  through  a  resistance 
Rj which  is  much  higher  than  the  transistor  resistance  This  decrease?  the  value 
of  ij  and  causes  it  to  be  essentially  in-pease  with  ?.  "he  actual  cirrent  passing 
through  the  winding.,  l,  can  be  or.  the  order  of  100  times  "he  current  ij  controlling 
the  transistor,  i  e.,  the  evans’seor  aits  as  a  current  amplifier  Due  to  the.  oper¬ 
ating  characteristics  of  transistors  "he  current  i  will  be  proportional  to  and  in- 
phase  with  i^(and  therefore  in  pbaee  with  1}  and  will  be  independent  of  the  voltage 
drop  L  di/dt.  provided  E,  >  >.L  di/dt.  +  R„i’-  »  AE  ,  where  i?  »,i<?  vol  tage  thres- 

fl  c  mm  i  tri<3  k 

hold  of  the  transistor.  A£  is  typically  1  \o  2  volts.  3nd  gen-ra'ly  ’ess  than  5 
volts  It  is  pointed  out  here  that  the.  power  loss  due  to  i^  R.^  will  be  relatively 
low  because  1^  is  on  the  order  of  i/.’GO 
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Since  the  driving  voltage  or.  EM  (T)  an d  Ed  (T)  -a  propor:  o -■«  i  *  :  t  h-  sn^t: 
velocity,  one  can  write  by  letting  me  proportionality  constant  be  y 

de 

_  x 

E  =  v  r 


-  -  yne  _  sin  fit 
xt 


Eq.  (43)  has  been  used  in  deriving  Eq.  ,’3C).  Note,  t  oat  E  is  a  lie  r  r.a  t  <=  1  ,■ 
applied  to  EM  ^  in  cne  interval  0  <  Ot  <  tt,  ace  to  EM  (5)  in  the  interval 
tt  <  nt  <  2tt.  The  electric  currents  in  the  windings  are  agai  i  alter r.a cel> 

E  %  dex  r 

1  s  “r  =  i  f  =  -  “sr 


in  the  respective  intervals.  R  is  the  external  rest stance  in  series  wi-h  cne 
winding,  or  in  the  case  of  the  tran-ntor  sir&uit,  an  equivalent  constant.  The 
combined  time-dependent  msg-’fe.. forte,  Ed  (T)  arc  F.d  (5)  by  virtue 

of  Eqs.  (15),  (16)  and  (51),  tan  be  expressed  a 


<*«>« 


5 _  N  y-,2 _ x  > 

-6  v0.626  i  R'  dt  dt 

10  g 


Employing  an  identical  arrangement  Tt.  the  y  dire  tt  ton  •  see  fig-  -  Ed  (F)  » r.c 
EM  (4)  would  similarly  exert  5  time  dependent  nag  e:  ic  f  orce  give:,  by 


<v« 


72  x  I06  °*626  *g  F 


v.  2  'Jj 

?'  dt 


Note  that  these  forces,  while  acting  to  oppose  the  velocity  of  tn-  snat: 
center,  are  proportional  to  tne  square  of  its  magnitude.  Therefore  tie  sta 
bility  analysis  presented  in  Section  II .  in.  which  various  force  compor.e-ts 
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are  rpnn-freH  iir.cirl”  prcpcrticr.cl  tc  either  un  uia^jiac emen c  or  the 

velocity  of  the  shaft  center,  can  not  be  applied  directly.  In  Appendix  IV,  '♦ 
a  virtual  damping  criterion  is  derived  showing  that  if  magnetic  forces  such 
as  those  given  in  Eqs .  (52)  and  (55)  are  employed  to  control  the  dynamic  or¬ 
bit  of  a  shaft-bearing  system,  which  would  be  otherwise  unstable,  an  orbit  of 
finite  amplitude  would  develop.  Thus,  the  type  of  magnetic  forces  discussed 
in  this  section  does  not  stabilize  the  shaft  in  the  absolute  sense,  but  in¬ 
stead,  it  limits  the  amplitude  of  the  shaft  orbit.  Equations  (52)  and  (53) 
can  be  rewritten  in  a  dimensionless  form  consistent  with  common  conventions 
employed  in  bearing  analysis,  e.g.  in  Section  II  and  Appendix  IV. 


<54) 


(55) 


Let  us  consider  for  example  the  case  of  an  unloaded  journal  bearing.  This 
bearing  is  inherently  unstable.  However,  in  the  presence  of  stabilizing 
forces  of  the  type  considered  here,  the  whirl  motion  is  limited  to  an  orbit 
of  finite  magnitude.  As  treated  in  detail  in  Appendix  IV,  the  dominant  com¬ 
ponent  of  such  an  orbit  as  estimated  by  a  one-step  iterative  Fourier  analysis 
would  be  a  circle  with  the  dimensionless  radius 


3tt  .^1  ~  V//l  jj_  3tt  3l  ~  V//l 

8  Vi  +  V  "  8  1  f2p 
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where  (U,yj,  v y  j>  )  arc  dyiamically  pert  urbea  beari-ig  forces  .  iee 

Appendices  I  and  IV)  at  the  fiequsnty  rat:. 


f 


which  is  approximately  0.5  regardless  of  the  precise  value  of  m.  Thus.  : f  the 
orbit  Size  is  to  be  limited  to  a  designated  amplitude.,  tne  required  value  of 
3  must  be  sufficiently  large. 


For  a  loaded  journal  bearing,  /  0,  the  analysis  of  Appendix  IV  would  have 
to  be  extended  to  allow  for  the  lack  of  symmetry  between  x  and  y  directions. 
For  a  heuristic  first  approximation.,  one  may  replace  Eqs.  (54)  by  their  dom- 
inant  Fourier  components  evaluated  at  an  acceptable  upper  bound  of  their 
amplitudes,  say  e^«  In  other  words,  one  may  make  the  approximation 


i  V  0 

xt  m 

mxx 

m 

(W  .  ) 

0  IV 

yt  m 
—  — 

myy 

—  mm 

where 


exp( if t) 


(56) 


( 5  ' ) 


and  thin  uoe  the  linear  stability  analysis  according  to  Section  11.  For  _:r- 
stance,  with  »  0.1,  an  unstable  system  (according  to  the  linear  stability 
analysis)  would  have  an  orbit  with  its  dominant  Fourier  component  larger  than 
0.1,  Conversely,  a  stable  system  would  have  an  orbit  with  ita  dominant. 
Fourier  component  smaller  than  0.1. 


If  one  strives  for  a  means  to  reduce  the  orbit  size  to  an  absolute  zero,  the 
magnetic  damping  must  be  not  only  phased  to  oppose  the  instantaneous  velocity, 
but  should  also  be  made  linearly  proportional  to  its  magnitude.  Two  ways  of 
achieving  linear  magnetic  dynamic  damping  appear  possible.  These  are  by  means 
of  a  square  root  function  generator  and  by  means  of  d-c  current  bias.  Each  has 
serious  shortcomings.  In  what  follows  below  their  operating  princ  iples  will  be 
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explained  and  shortcomings  pointed  ojt,  so  that  one  car;.,  in  the  future  r  ;ii- 
sider  them  as  alternatives  ins  more  rigorous  trade-off  study. 


Square-Root  Function  Generator 

A  square- root  function  generator  (SRFG,  is  a  device  which  when  fed  by 
an  input  f  (t)  wo  uld  produce  an  output  .  If  we  connect  in  aeries 

one  SRfG  before  each  electromagnet ,  then  it  is  quite  obvious  that  the 
force  produced  will  be  proportional  to  e  .  This  will  enable  us  to 
use  the  linear  stability  theory. 

One  such  SRFG  ve  have  found  la  the  General  Electric  Square-Root  Con¬ 
vertor  (GEMAC  System  HBK  8058fi>.  According  to  the  G.E.  Brochure,  it 
has  the  following  characteristics: 


Weight 


5  pound  net 
3"  x  6"  x  15” 


Frequency  Response 


Down  3db  at  3  Hz 


The  frequency  response  is  definitely  unsatisfactory  because  we  are  in¬ 
terested  in  the  range  of  100  Hz  or  higher. 


D-C  Current  Bias 

Apply  a  d-c  bias  electric  current  m  the  coll  of  each  electromagnet. 

\ 

This  bias  current  should  be  at  least  one  order-of  magnitude  larger  than 
the  current  driven  by  the  amplifier.  |Let  and  1^  denote  respectively 

the  d-c  bias  electric  current  and  the 'current,  driven  by  the  amplifier 

i 

(I-  >  >  I1).  Now  the  force  produced  by  the  electromagnet  is  proportional 

coUq+V2.  j 

2  2 

Electromagnet  Force  —  ( IQ  +  1^  e*  7.^  +  2T(.T^.  Note  that  we  can  neglect 

the  Ij  term  with  the  use  of  the.  d-c  bias  current.  The  force  produces  by 
the  electromagnet  is  therefore  proportional  to  I  and  hence  e  .  Again., 
one  can  now  use  the  linear  stability  theory. 
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However,  this  bias  >;-;trert  giv<  -  r.s-  ■  :>  a  steady  magnetic  force  which  in 
turn  results  in  a  negative  stiffness.  Inis  negative  stiffness  should  oe 
appreciable  because  ?  ,T,  .  A.lsc  the  application  of  tfus  bias  current 

causes  a  large  increase  in  power  loss  which  will  be  discussed  in  the  sec¬ 
tion  entitled,  Power  Loss  and  Weight  Penalties. 

In  summary,  we  do  not  recommend  the  use  of  either  the  SRFG  or  the  d-c  Bias  Current 
just  for  the  purpose  of  making  the  '.inear  stability  theory  applicable  Due  to  the 
lack  of  knowledge  about  the  influence  of  non-linear  dynamic  damping  it  is  difficult 
to  predict  whether  a  linear  or  non-linear  dynamic  damping  is  more  effective  in  sup¬ 
pressing  whirl  instability.  But,  qualitatively  it  is  beyond  doubt  that  the  non¬ 
linear  dynamic  damping  produced  by  the  proposed  active  device  will  be  effective  in 
suppressing  instabilities.  In  the  following  example,  some  preliminary  stability 
results  will  be  obtained  by  using  the  linear  stability  theory  and  Eq .  (57) 

EXAMPLE:  Loaded  Plain  Journal  Bearing 

/ 

Input 

D  ■  1  inch 

C  ■  0.001  inch 

u  »  0.27  x  10  ®  lb.sec/in^ 

A  -  0.35 

e  -  0.2 

xo 

L/D  *=  1 

From  the  definition  of  A  we  calculate  the  journal  speed  u>  «  1270  rad/sec.  If  we 
further  specify  the  following: 

Additional  Incut 
N  ■  30  turns 

i  *  1  amp 

R  *=  10  ohm  (including  external  resistance) 
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»  11,8  volt 

A  -  0.2  inZ 

S 

£  =  0.01  Inch 

8 


In  order  to  have  E ^  *  11.8  volts,  and  suppose  that  the  tolerable  la  0.118, 
then  it  is  required  that  the  amplifier  can  yield  y  fl  C  *  100  volts.  With 
the  above  input  data,  the  magnetic  dynamic  damping  can  be  readily  calculated 
from  (57). 


V  -  V 
irocx  myy 


0,0139 


The  values  of  V  and  V  can  be  easily  increased  by  increasing  the  ampli- 
mxx  myy  J  J 

fication  ratio  of  the  amplifier,  y„ 


Before  we  go  any  further,  let  us  check  if  inequality  (48)  is  satisfied.  First 
of  all,  ft  is  typically  one  half  (or  smaller)  the  Journal  speed,  ft  »  635  rad/sec. 
From  Eq.  (23),  L  ■  2.87  x  10*^  henry 

a.  h?..»  2.87 .  iP;4  .  0  0l82 

which  is  indeed  much  smaller  than  tt  required  by  inequality  (48), 

The  linear  stability  results  for  different  values  of  magnetic  dynamic  damping 
are  tabulated  as  follows: 

y  (.  V  )  Dimensionless 

mxx  ^ *  bw _ Threshold  Speed  (Eq.  12a) 


0 

0.286 

0.005 

0,375 

0.010 

0.575 

0.015 

0.891 

0.020 

1.373 

At  V  ■  0.0139,  the  threshold  speed  is  0,  83,  which  is  three  times  the 

value  at  V  ■  0. 
mxx 
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POWER  LOSS  AND  WEIGHT  PENALTIES 

In  the  previous  sections  we  have  shown  that  both  the  device  producing  unidi¬ 
rectional  magnetic  force  and  the  active  device  are  capable  of  making  a  journal 
bearing  more  stable.  It  is  then  m  order  to  estimate  what  penalties  are 
associated  with  the  use  of  either  of  the  devices  One  major  penalty  in  of 
course  the  power  loss  associated  with  tne  electromagnetic  sysiem  The  other 

major  item  is  the  weight  penalty  if  the  eae  bearing  is  to  be  used  in  a  space 

unit  or  any  other  circumstances  where  weight  is  an  important  consideration. 

It  is  illustrative  to  investigate  an  actual  machine  utilizing  gss  bearing 
suspension  and  compare  the  power  loss  with  the  machine  power  rating  and  the 

additional  weight  with  the  total  weight  of  the  machine. 

Tilting-pad  gas  bearings  supporting  a  9  kw  turboalternator  were  designed 
by  MTI  for  NASA  (under  subcontract  to  Pratt- Whitney  Aircraft)  The  bearings 
had  been  successfully  operated  for  more  than  1200  hours  and  they  are  still  in 
operation.  The  bearing  dimensions  and  generator  characteristic  data  are 
listed  below  for  easy  reference. 

9  kw  ,  4  pole .  400  Hz 
Shaft  diameter  -•  3-1/2  inch 

Shaft  speed  **  12  000  rpm 

.  Stator-rotor  sir  gap  *•  20  ml  1 

Rotor  mass  56  lb 

1.  Preliminary  Design  -  Unidirectional  Device 

Suppose  that  we  want  to  operate  the  unit  described  above  in  a  zero-g  environment 
with  two  plain  journal  gas  bearings  whose  D  •  3-1/2  Inch,  C  *  0  001  inch  and 
L/D  “  1.  One  possible  way  to  achieve  a  stable  operation  is  utilizing  the  de¬ 
vice  producing  a  unidirectional  magnetic  force.  In  order  to  obtain  a  critical 
mass  of  28  pounds  which  is  half  of  the  rotor  mass  (because  there  ere  two 
bearings),  It  is  found  that  the  bearings  have  to  be  loaded  magnetically  to 
an  eccentricity  of  *  *  0,41,  or  F  *-  82  lb„  Let  the  flux  density 

m  XO  TT1 

be  40,000  line/in  ,  Then,  the  area  of  the  air  gap  needed  is 
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r  x  72  x  106 
.  m 

A»  * - 5 - 

B" 

Design  a  cross- sect! or.  of 
Assume  an  air  gap  of  /  ■ 

gap  i», 


82  x  72  x  106  .  .  2 

—  •  —.1  ■— .  '•  m  3.69  in 

(40, 000) ‘ 

,  say,  v  m  1-7/B  inch  and  d  =  2  inch  (see  Fig.  13) 
0.01  inch,  the  ampere- turn  required  for  the  air 


(Hi)  “  (0.313  B  X  )2  -  250  amp-turna 

gep  S 

Assume  50  amp- turns  to  take  into  account  the  drops  in  iron  and  leakage.  Thus 

to  *  250  +  50  *  300  amp- turns.  Choose  i  •  1  sap,  N  •  300  turns  and  use  #18 

wire  which  has  an  o.d.  of  0,04  inch  and  A.  -  1600  c.m.  The  cross-section  of 

K  2 

the  coll  winding  would  be  0.04  x  0,04  x  300  »  0.48  in  .  Allow  an  additional 

o 

(see  Fig.  13)  0.75  in  for  insulation.  Thus,  to  acconmodate  the  coil  winding 

2 

and  some  margin  we  design  a  cross-sectional  area  of  2  in  ,  ao  that 


Xy  K  2  inches,  «  1  inch 

2.  Preliminary  Design  -  Active  Device 

Design  the  system  to  be  able  to  yield  a  maximum  of  V  ■  0.1  at  *  0,118. 

Then,  the  peak  magnitude  of  magnetic  time- dependent  force  is 

(p  tt  D  1)  V  «, 
w  '  mxx  1 

«  (14,7  x  TT- x  3.5  x  3.5)  0.1  x  0.118 

*»  6. €7  pounds 

o 

Assume  B  «  20,000  llnes/in  .  Then,  the  area  of  the  air  gap  ia 

6.67  (72  x  106)  .  „  .  i 

A  «  - * - - — —  -  -  1.2  in 

8  2C  s000 

We  car.  make  the  croea-secticn  to  be,  say,  w  -  1.1  inch,  d  •  l.l  inch.  Again 
assume  an  air  gap  of  0.01  inch.  The  empere-turn  required  is  125.  Allow  an 
additional  25  amp-turna  for  leakage,  etc.  We  need  a  peak  total  of  150  amp-turns. 
Chooae  a  peak  current  of  1  amp  ar.d  uae  #18  wire.  H  *  150  turns. 


-32- 


Since  tne  ampere  turn  is  one  half  of  that  of  the  unidirectional 

device  the  cross-sectional  area  of  the  coil  and  insulation  is  also  about 

half  of  it.  (0.48  +  0.75)/2  -  C.62  in*. 

a.  Loaded  Bearing 

If  the  bearing  is  loaded  by  its  own  weight,  (28  lb/bearing)  it  will  reach  a 
static  equilibrium  position.  For  C  =  0.001  inch,  the  bearing  will  develop 
this  load  capacity  at  exo  «=  0.15.  The  ctitical  mass  can  be  calculated  to 

be  three  pounds.  Since  the  critical  mass  is  smaller  than  one  half  of  the 

•  <- 

rotor  mass,  the  bearing  would  be  unstable.  In  order  to  uHfct  the  rotor 

bearing  system  stable,  the  critical  mass  should  be  no  less  than  28  pounds. 

This  can  be  achieved  by  using  the  active  device.  A  plot  of  the  critical 

mass  against  V  is  shown  in  Fig.  11.  It  is  seen  that  a  V  of 
mxx  mxx 

approximately  0.103  Is  needed  to  achieve  a  critical  mass  of  28  pounds. 

In  fact,  if  we  design  the  bearing  to  have  a  larger  clearance,  the  bearing 

would  operate  at  a  higher  but  have  smaller  dynamic  stiffness.  As  a 

result,  the  critical  mass  would  have  a  lower  value  when  V  »  V  »  0. 

mxx  myy 

However,  it  rises  rather  quickly  with  increasing  magnetic  damping  as  shown 
in  Fig.  11.  The  following  table  summarizes  the  magnetic  damping  needed  for 
various  C  to  achieve  a  critical  mass  of  28  pounds. 


C  [in] 

A 

c 

XO 

V 

mxx 

Critical  Mass  [lb] 

0.0010 

4.28 

0.15 

0.103 

28 

0.0015 

1.90 

0.24 

0.092 

28 

0.0020 

1.07 

0.39 

0.076 

28 

b,  Unloaded  Bearing 

If  the  bearing  is  unloaded  (e.g.,  m  a  zero-g  environment),  the  critical 

mess  is  again  plotted  against  V  x  in  Fig.  12.  Note  that  the  critical  mass 

approaches  zero  when  V  -  0.  This  Indicates  the  well-known  fact  that  an 

mxx 

unloaded  Journal  bearing  is  always  unstable.  When  *  0.1,  the  critical 

mass  becomes  as  high  as  252.6  pounds,  in  the  following  we  again  summarize 
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the  dynamic  magnetic  damping,  V  ,  needed  for  various  C  to  achieve  a  critical 
ma3s  of  28  pounds  8t  zero  eccentricity. 


C  [in] 

V 

mxx 

Critical  Maaa  [lb] 

0.0010 

0.109 

28 

0.0015 

0.096 

28 

0,0020 

0.077 

28 

It  ia  therefore  recommended  to  design  the  bearing  at  C  ■  0,00?  Inch,  if  the 

active  device  is  to  be  used.  To  operate  stably,  it  would  need  V  ■  V  • 

..  mxx  myy 

0.076  when  it  is  loaded,  V  »  V  ■  0.077  when  It  is  not  loaded.  For  sub- 

mxx  myy 

sequent  calculations,  0,08  will  be  used  for  the  value  of  V _  .  In  Section 

mxx 

IV- 7,  it  will  be  seen  chat  larger  clearance  would  result  in  lower  bearing 
frictional  loss  although  the  eccentricity  will  be  higher. 

3.  Eddy  Current  Loss 

The  magnetic  flux  while  passing  through  the  shaft  are  cut  by  the  rotation  of 
the  shaft  and  thus  generating  e.m. f.  This  e.m.f.  will  produce  eddy-current 
and  dissipate  in  the  form  of  heat.  In  the  following  we  will  make  an  estimate 
of  this  eddy-current  loss. 

Magnetic  flux  will  complete  its  path  from  one  air  gap  through  the  shaft,  the 
other  air  gap  and  back  into  the  magnet.  It  will  penetrate  the  entire  depth  of 
the  shaft  in  a  three-dimensional  fashion  which  will  be  affected  by  the  ratio  of 
the  magnet  dimensions  to  the  shaft  diameter,  the  variation  in  permeability  due 
to  change  in  3,  and  the  preaence  of  eddy-current  in  the  shaft. 

While  an  actual  analysis  for  the  flux  distribution  is  rather  complicated,  we 
assume  for  a  first  approximation,  that  the  flux  is  flowing  uniformly  in  a 
channel  inside  the  shaft.  The  channel  penetrates  to  depth  d  and  width  w  in 
the  shaft  (i.e.,  the  same  dimensions  aa  the  outer  core)  as  shown  in  Fig.  13. 

If  the  rotation  of  the  shaft  la  aa  shown  in  Fig.  13,  and  the  flux  direction  is 
to  the  right  aa  shown,  then  the  e.m.f.  generated  ia  radial  in  the  central 
region  and  axial  near  the  gap  according  to 

d  Ein  *  'v  *  B  df)  10  0  (58) 
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A  current  (eddy-current)  Is  produced  because  of  the  induced  e.m.f.  and  is 
in  the  same  direction  as  E  .  A  force  proportional  to  the  product  of 
E  -rd  I.UUCUI.,  accs  m  a  cirection  to  oppose  the  shaft  motion.  Hence, 

there  is  a  drag,  and  eddy-current  loss  must  be  supplied  mechanically  by  the 
shaft . 


Let  us  focus  our  attention  in  the  central  region  where  the  flux  is  toward 

the  right  and  Is  in  the  radial  direction.  The  magnitude  of  the  velocity 

is  v  *  cur  and  At  ■  dr  where  co  is  the  rotating  speed  of  the  sheft  and  r  the 

radial  coordinate.  Using  the  above  and  integrating  (58)  we  have 

D 

2 

,-8 


in 


->»  5 


r  dr  x  10 


-  B  05  (Dd  -  d2)  x  10"8 
Substituting  values  from  the  preliminary  design 

B  •  40,000  lines/in^ 
cu  -  12000  rpm  -  1255  rad/sec 
D  ■  3.5  inches 
d  ■  2  inch 

_E  "  2  X  40»000  x  1255  (7  —  4)  x  10"8 

■  0.756  volts. 

In  order  to  calculate  the  eddy-current  loss  due  to  this  induced  voltage, 
it  is  necessary  to  estimate  the  electrical  resistance  of  the  current  path. 
Assume  that  current  is  flowing  over  an  area  of  cross  section  equal  to  the 
length  and  width  of  the  magnet.  Thus, 

Area  of  cross  section  ■  (£  +  d)  w 

w 

Length  of  current  path  **  d 
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The  electrical  resistance  of  the  induction  path  *  23.6  x  10 


(}  * 
'  w 


-  23.6  x  10‘ft  x 


,-6 


(2  4-  2)1.88 


*  6.27  x  iO  ohm  if  the  shaft  is  made  of  super* 


malloy  which  has  a  resistivity  of  23.6  x  10  ohm  per  inch  length  per  square 
inch  cross  section.  (See  Page  5-180  of  Ref.  14).  Since  electrical  current 
will  return  through  some  path  in  the  shaft,  we  assume  that  the  total  resistance 
is  twice  the  value  calculated  above 

2 


E 


R  --  12-5  x  10  ohm.  Eddy-current  power  loss 
watts  *  45 -6  kw . 


in 

R 


0 . 756 


12.5  x  10 


-6 


-  45-6  x  1CT 


Thus,  the  estimated  eddy-current  power  loss  is  about  5  times  greater  than  the 
power  rating  of  the  turboalternator .  Whatever  the  assumptions  chat  are  made 
regarding  the  flux  distribution  and  current  path  etc.,  it  is  obvious  that  the 

a 

eddy-current ‘ losses  will  be  prohibitive  at  the  speed  and  flux  density  specified. 

If  we  use  che  active  device  instead,  the  eddy-current  loss  is  then  estimated  to 
be  i  kw  which  is  still  very  high. 


A  common  practice  in  electrical  machinery  design  to  reduce  eddy-current  loss  is 
by  lamination,  but  it  appears  very  difficult  to  laminate  a  shaft.  We  there¬ 
fore  propose  to  mount  a  laminated  disk  onto  the  shaft  as  shown  in  Figs  14 
and  15.  Note  that  the  electromagnets  in  Figs.  14  and  15  have  been  oriented 
ar  such  a  way  that  the  induced  e.m.f.  is  always  perpendicular  to  the  lamina¬ 
tions-  The  laminated  disk  and  the  shaft  are  separated  by  a  ring  made  of  non¬ 
magnetic,  nonconducting  material  such  as  ceramic.  It  is  seen  that  the  induced 
*  :n  f  is  now  only  in  the  regions  near  the  gaps-  Suppose  that  there  are  q 
laminations  The  induced  e.m.f.  and  electrical  resistance  in  each  lamina¬ 
tion  are  E /q  and  qR,  respectively.  The  total  eddy-current  loss  is  clearly 
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1 


2 


n 


(hJly 

<1  R 


£. 

M 


Thus,  the  eddy-current  loss  being  proportional  to  the  inverse  square  of  q5 
can  be  very  effectively  reduced  to  a  tolerable  level  by  the  addition  of  a 
laminated  disk.  A  0„02  inch  thickners  lamination  is  quite  common.  It  will 
be  seen  later  that  the  dii>k  thickness  is  designed  to  be  four  inches  for  the 
unidirectional  device  and  one  inch  for  the  active  device. 


Thus, 


q  *  cT02  "  laminations  for  the  unidirectional  device 


q  *  "o“02  *  50  laminations  for  the  active  device 


.  45,600  watts  45.600 

Eddy-current  loss  »  ■* . r  - - -  -  " 

q  200 


1.1  watts  for  the  unidirectional  device 


Eddy-current  losa 


4000  watt a 
2 


4000 

"so2 


-  1.6  watts  for  the  active  device 


4.  Proposed  Design 

Since  it  is  always  desirable  to  have  the  magnetic  forde  concentrating  at  a 
particular  radial  direction  rather  than  over  a  wide  area,  a  bettor  magnet  de¬ 
sign  is  shown  in  Fig.  16  in  which  the  spacing  between  the  pole  pieces  is  re¬ 
duced.  Data  of  the  improved  designs  for  the  magnet  and  the  laminated  disk  are 
tabulated  in  the  following.  The  subsequent  estimates  on  power  losses  and 
weight  penalty  are  based  on  this  configuration. 
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Unidirectional 

Device 


Active  Device 


laminated  Disk  Thickness 

O.D. 

I.D. 

Magnet  face  Width,  d 
Magnet  Face  Depth,  w 


Coil  Cross-sectional  Area 

*w 


Number  of  Magnets  Needed 


4  inches 
6  inches 
4  inches 
0.92  inch 
4  inches 

1.25  in2 
1.25  inch 

1  inch 

1 


5.  Hysteresis  Loss 

The  hysteresis  loss  is  given  by  (Ref.  6) 

\  ’  \  c  \  1,4 


1  inch 
6  inches 
4  inches 
1.1  inch 
1.1  inch 

0.62  in2 
1  inch 

0.62  inch 

4 


2 

where  f  us  the  frequency  in  Hz,  B  is  the  maximum  flux  density  in  lines/in  , 

ID  X3 

and  K.  ;.a  a  material  constant  (which  for  the  units  stated  above  is  4  x  10 
n 

for  supermalloy,  Ref.  13,  page  513). 


In  the  device  producing  a  unidirectional  magnetic  force,  the  flux  density, 

2 

according  to  the  preliminary  design  is  40,000  lines/ln  .  The  flux  fluxuation 

2 

ts  due  to  shaft  rotation,  so  that  B  ■  40,000  lines/in  and  f  *  2v0  cps.  In 

m  2 

Che  active  device,  B  is  20,000  lines/in  and  the  frequency  f  is  100  cps.  if  we 
tn 

aasume  a  ha If- frequency  whirl  motion.  Thus,  the  hysteresis  losses  are,  using 
supermalloy  as  the  magnetic  material, 

F,  =  4  x  10'13  x  100  x  20.0001*6  -  0.0003  watts/lb 

n 

for  the  active  device,  and 

- 13  1  6 

Ph  »  4  x  10  x  200  x  40,000  *  -  0.0018  watts/lb 

for  the  unidirectional  device. 
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The  elements  having  magnetic  flux  fluctuation  are  the  laminated  disk  and  magnet 
cores  in  the  active  device,  the  laminated  disk  in  the  unidirect ional  devfrp  A« 
shown  in  Section  IV-9,  the  weight  of  these  elements  is  4,98  +  4  x  2.92  -  16 .7  lb, 
for  the  sf five  device,  and  19-9  lb.  for  the  unidirectional  device.  Thus ,  the 
hysteresis  ios3  is 

0.0003  x  16.7  ■  0.  i  watts  for  the  active  device  and 
0.C018  x  19.9  --  0.036  watts  for  the  unidirectional  device 


6.  Copper  Loss 


a.  Unidirectional  Device 

Mean  length  of  one  turn  ■*  (w  +  +  d  +  /^)2  ”  (4  +  1  +  0.92  +  1)2 

-  13.8  in. 

Total  length  **  300  x  13.8  *  4130  inches  »  344  feet 
Area  of  #18  copper  wire  «  1600  c.m. 


R 


10.7  x  3.44 
1600 


2.3  ohm 


i  »  current  »  1  amp 

2 

Power  loss  -  i  R  -  2.3  watts 

2  3 

For  each  bearing,  the  copper  less  is  -j-  ■  1.15  wstts/bearing. 


b.  Active  Device 

Mean  length  of  one  turn  *  (w  +  +  d  +  1^)2  ■  (1.1+0.62+1.1 

+  0.62)2  *  6.88  in 

Total  length  »  150  x  6.88  ■  1030  inches  ■  86  feet 
Area  of  #18  copper  wire  •*  1600  c.m. 


R  « 


10.7  x  86 
1600 


C.575  ohm 


As  indicated  before,  an  external  resistance  is  often  necessary  in  order  to  make 
the  circuit  resistive  rather  than  inductive.  The  inductance  is  by  Eq.  (23). 


L  «  3.19  x  10 


3.19  x  10‘8  -y—-*  KT±  *  4.35  x  !0'2  henry 

2  x  0.01 


Assume  that  the  whirl  speed  is  half  the  shaft  speed,  -:hen 


W  =  —  1225  =  628  rad /sec 


QL  =  628  x  4.35  x  10 


27.3  ohm 


If  we  add  an  external  resistance  to  make  a  total  resistance  of  “  140  ohm,  then 
QL/R  *  0.195  which  i9  one  order  of  magnitude  smaller  than  u  Wi'ts  a  peak  current 
of  1  amp,  the  current  in  EM  is,  from  Eqs  .  (46)  and  149).  i.  s‘  i  sin  fit. 


Thus, 


r 

copper  loss  in  EM  \  t^R  d  (Cit)s-r^  \  t2  sin2  fit  d  (fit) 

c  2T»  max 


^max^e  140 


■  35  watt 


Since  there  are  four  electromagnets  in  the  active  device,  total  copper  loss  * 
4  x  35  *  140  watts-  And  because  there  are  two  bearings,  t.hu  copper  loss  per 
bearing  is  140/2  ■  70  watts/bearing. 


As  Indicated  on  Page  23,  one  can  use  a  transistor  circuit  ins*,  cad  of  adding  external 
resistance,  to  achieve  the  same  purpose.  From  the  calcul-.r  i  r.r.s  ui  pages  .39  and  40, 


we  have 


L  -  4.35  x  10  henry 

R  *  Resistance  of  coil  winding 

*  0.575  ohm 

i  =*  I  sin  fit 
max 

1=1  amp. 
max  r 


ThuSj  we  calculate 


J\ _ gear  ini'  Film  Loss 

The  loss  in  the  bearing  film  is  predominantly  frictional  loss  due  to  shear. 
ih«  «hear  stress  is  equal  to  the  product  of  viscosity  and  velocity  gradient. 


■sm). 

c 


tne  effect  of  eccentricity  has  been  neglected. 


To  obtain  power  loss,  we  should  multiply  the  shear  stress  by  ar.  area  (rrDL). 
and  a  velocity  ujD/2.  Thus, 


Bearing  film  loss 


u  (ttDL)  tip 


With 


M  *•  2.7  x  10  9  lb-sec/in^ 


OJ 

D 


C 


1255  rad/sec 
L  *  3.5  inches 

0.001  inch  for  unidirectional  device 

0.002  inch  for  active  device 


we  calculate 


Bearing  film  loss  <*  < 


56.5  watts/bearing  for  unidirectional  device 


28.3  watta/bearing  for  active  device 


8.  Total  Power  Loss 

The  total  power  loss  is  the  sum  of  eddy-current  loas,  hysteresis  loss,  copper 
lo-is  and  bearing  film  loss.  Thus, 

Total  Power  Los3  *  1.1  +  0.036  +  2.3  +  2  x  56.5 

*  116.4  watts  for  unidirectional  device 

Total  Power  Loss  *  1.6  +  0.005  +  140  +  2  x  28.3 

*  198.2  watts  for  active  device  with  external  resl8tance 
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Total  Power  Loss  =  1.6  +  0.005  +  41.2  +  2  x  28.3 

«■  99.4  watts  for  active  device  with  transistor  circuits 

The  above  power  losses  are  to  be  compared  with  the  loss  in  the  tilting  pad  gaq 
bearings  currently  in  use;  it  was  reported  to  be  180  watts  for  the  two  bearings 


9.  Weight  Penalties 

a.  Unidirectional  Device 

The  laminated  disk  has  a  thickness  of  4  inches,  and  outside  and  inside 
diameters  of  6  in.  and  4  in.  respectively  (See  Section  IV-4) . 


Volume  of  Laminated  Disk  *  4  x  £  (62-  42)  ■  62.8  in3 
Weight  of  Laminated  Disk  ■  Density  x  Volume*  0.317  x  62.8 

■  19.9  lb.  (Supermalloy) 


3 

The  Density  of  Supermalloy  is  0.317  lL/in 

Volume  of  Magnet  Core  ■  (v  x  d)  (jly  +  d  +  +  d)2 


Weight  of  Magnet  Core 


(4  x  0.92)  (1.25  +  0.92  +  1  +  0.92)2 
30.3  in2 

Density  x  Volume  *  0.317  x  30.3  *  9.6 


lb. 


Volume  of  Coil  Winding 


(Supermalloy) 

(I  x  V)  (d  +  V  +  w  +  £' )2 
www  w  3 

(1.25  x  1)  (0.92  +  1  +  4  +  1)2  -  17.3  in 


From  Section  IV- 1,  the  cross  sectional  areas  of  the  copper  and  the  insula- 

2  2 

tion  in  the  winding  are  respectively  0.48  in  and  0.75  in  .  Therefore, 
the  volumetric  percentage  of  copper  in  the  winding  is  0.48/(0.48  +  0.75) 

*  39%,  and  the  remaining  61%  is  insulation.  The  density  of  copper  is 

3  3 

0.324  lb/in  .  Assume  the  density  of  insulation  to  be  0.065  lb/in  (20% 

of  copper  density).  Then,  the  density  of  the  coil  winding  is  (0.324  x 

0.39  +  0.065  x  0.61)  -  0.166  lb/in3. 


Weight  of  Coil  Winding 


■  Density  x  Volume  *  0.166  x  17.3  ■  2.87  lb. 
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Since  only  a-c  power  is  available  from  the  turboalternator  and  d-c  power 
is  needed  for  the  electromagnet,  a  rectifier  is  required  whose  weight  is 
estimated  at  0.5  lb. 

Total  Weight  Penalty  -  19.9  +  9.6  +  2.87  +  0.5  *  32.9  lb. 

b.  Active  Device 

The  laminated  disk  has  a  thickness  of  1  inch,  and  outside  and  inside 
diameters  of  6  in.  and  4  in.  respectively  (See  Section  IV- 4) . 

Volume  of  Laminated  Disk  “  1  x  j  (6^  -  4^) 

■  15.7  in^ 

Weight  of  Laminated  Disk  ■  Density  x  Volume  ■  0.317  x  15.7  *  4.98  lb. 

(Supermalloy) 

.j.  V.oliime  of  One  Magnet  Core  ■«  (w  x  d)  ( i ,  +  d  +  +  d)2 

■  (l.lx  1.1)  (U  1.1  +  0.62  +1.1)2 

-  9.2  in3 

Weight. of  One  Magnet  Core  •  Density  x  Volume  *  0.317  x  9.2 

■  2.92  lb.  (Supermalloy) 

Volume  of  One  Coil  Winding  -  (i  x  V)  (d  x  l'  +  w'+  l')2 

w  w  w  w 

-  (1  x  0.62)  (1.1  +  0.62  +  1.1  +  0.62)2 

-  4.26  in3 

3 

The  density  of  coil  winding  is  as  shown  in  Section  IV-9-a,  0.166  lb/in  . 
Weight  of  One  Coil  Winding  ■  Density  x  Volume  »  0.166  x  4.26  ■  0.708  lb 
Total  Weight  of  Disk,  and  Four  Sets  of  Magnet  and  Coil  Winding 

-  4.98  +  4  (2.92  +  0.708)  *  19.5  lb. 
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Allow  three  pounds  for  the  amplifier-transformer  unit  and  0.5  lb,  for  the 
proDe,  ditterentiator  and  signal  splitter  (See  Fig,  8);  and  there  are  two 
subsystems  in  the  active  device. 

Total  Weight  Penalty  *  19.5  +  2  (3  +  0.5)  *  26.5  lbs. 
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SUMMARY 


The  major  results  of  this  feasibility  st-ndy  are  summarized  as  follows; 

1.  Two  devices  appear  to  be  feasible  to  Improve  the  stability 
of  rotor-bearing  systems. 

A.  Unidirectional  Device--producing  a  unidirectional 
magnetic  force  to  load  the  bearing. 

B.  Active  Device--producing  a  controlled  electromagnetic 
force  always  opposing  the  motion  of  the  shaft. 

2.  Plain  journal  bearings  operating  in  a  zero-g  environment  dan  be 
made  stable  by  using  either  the  unidirectional  device  or  the 
active  device;  the  stability  margin  achieved  by  the  unidirectional 
device  is  at  best  equal  to  that  of  a  gravity- loaded  bearing  if 

A  is  large  or  moderately  large  (A  J>  0.3  roughly). 

3.  The  active  device  is  very  effective  in  suppressing  whirl 

Instabilities.  A  nominal  magnetic  damping  produced  by  the 

active  device  can  increase  the  critical  mass  several-fold. 

Preliminary  calculation  shows  that  for  A  *  1.0,  the  critical  mass 

is  increased  100  times  with  a  V  »  V  ■  0.1  which  can  be  easily 

mxx  myy  7 

achieved.  The  effect  becomes  even  more  profound  for  smaller  A, 

but  less  for  larger  A. 

4.  The  magnetic  damping  produced  by  the  active  device  is  not  linear 
with  the  rotor  displacement.  A  quasl-linear  approach  is  used  in 
the  critical  mass  calculation. 

3.  There  are  two  possible  ways  to  make  the  magnetic  damping  linear. 

One  is  to  use  square  root  function  generators,  and  the  other  is  to 
apply  a  d-c  current  bias.  Both  methods  have  their  shortcomings  as 
stated  in  the  text.  We  do  not  recommend  to  use  either  of  the  two 
until  it  is  proven  that  the  linear  magnetic  damping  is  auperior  to 
the  non-linear  one.. 
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6  A  9  kw  turboalternator  supported  by  two  tilting  pad  gas  bearincs 
was  used  as  a  reference  to  assess  the  power  loss  and  weight  penal¬ 
ties  associated  with  the  use  of  either  devicp ,  They  are  listed  in 
the  following. 

Plain  Journal  Bearings  with  Unidirectional  Device 

Power  Loss  *  116.4  watts 

Weight  Penalty  =  32.9  lbs. 

Plain  Journal  Bearings  with  Active  Device 

Power  Loss  *  99.4  watts 

Weight  Penalty  »  26.5  lbs. 

Comparing  with  the  rotor  weight  of  56  lbs.,  the  weight  penalty  of 
using  either  device  is  appreciable.  The  above  power  losses  compare 
favorably  with  the  existing  bearing  system  as  the  frictional  loss 
in  the  two  gas  bearings  of  tilting-pad  design  was  reported  to  be  180 
watts . 
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RECOrififcaUAilUHS 


It  Is  recommended  to  design  and  fabricate  a  rotor-bearing  unit  employing 
plain  Journal  gas  bearings  and  the  active  electromagnetic  device.  The 
rotor-bearing  system  should  be  designed  to  be  in  a  dynamically  unstable 
condition  without  the  aid  of  the  active  device.  By  switching  the  active 
device  on  and  off,  one  would  be  able  to  demonstrate  its  effectiveness  in 
suppressing  the  whirl  instability  of  the  rotor. 
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Fig.  3 


Magnetic  Stiffness  and  Damping  versus  Q  L  /R 
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Critical  Mass  versus  Dynamic  Magnetic  Damping 
under  Gravitational  Loading 


ASSUMED  FLUX  CHANNEL 


L-  d--a 


Unidirectional  Device  with  a  Laainated  Disk 


Fig.  16  Configuration  of  Electromagnet  in  Association 
with  a  Laminated  Disk 


/ 


In  Appendix  IV,  where  harmonica  of  the  whirl  frequency  ere 
following  nomenc leture  ia  uaed: 


considered. 


the 


U  +  IV 
xxn  iVxx„ 


n  being  the  harmonic  number,  n  «  1,  2,  3f 
and  n  *  1, 


In  particular,  for  «  .  0 

xo 
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APPENDIX  II;  Differentiator  and  Signal  Splltcer 


in  an  active  device,  differentiators  and  signal  splitters  are  required  as  in¬ 
dicated  in  the  text-  Their  operating  principles  and  circuit  diagrams  will  be 
illustrated  here. 

A  standard  differentiator  circuit  shown  below  can  be  assembled  to  perform  the 
differentiation  operation  upon  the  input  voltage  which  in  this  particular  ap¬ 
plication  is  the  output  of  the  displacement  probe.  In  the  sketch,  A  is  an 


C 

■U ■ 


AAAAr 


J input 


*> 


output 


T 


Differentiator  Circuit 


operational  amplifier,  C  a  capacitor  and  R^,  a  resistor.  The  resistor  and  the 
operational  amplifier  form  a  feedback  network  to  regulate  the  gain.  The  capac¬ 
itor  and  this  feedback  network  would  then  perform  the  differentiation  operation. 
A  more  detailed  explanation  on  the  differentiator  circuit  can  be  found  in  [8]. 


-69 


A  signal  splitter  consists  oi  two  diodes  and  as  shown  in  the  sketch. 


Signal  Splitter 

The  two  diodea  are  connected  so  that  a  positive- going  signal  will  be  routed  to 
E^t  and  a  negative -going  aignal  will  be  routed  to  Eout  The  two  output  term¬ 
inals  are  to  be  connected  to  the  two  coil  windings  of,  say,  EM  and  EM 
of  Fig.  8. 
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APPENDIX  III:  Computer  Program  -  PN  424  -  Influences  of  Magnetic  Forces  (Uni¬ 
directional  or  by  an  Active  Device)  on  Stability  of  Plain  Journal 
Bearings 

INPUT 

Card  1  (80H)  Title  card 


Card  2  (8EI0.3) 


Card  3  (8E10.3) 


Card  4  (8E10.3) 


Card  5  (415) 


This  card  contains  the  following  five  values: 

X  ■  magnetic  air  gap  [in] 

8 

C  *  radial  bearing  clearance  [in'! 

2 

u  =  viscosity  [lb-sec/in  ] 

p  *  ambient  pressure  [lb/in  ] 

d 

D  *  diameter  of  shaft  [in] 

This  card  contains  the  following  three  dimension¬ 
less  quantities: 

f  ■  initial  value  of  frequency  ratio  f 

*  final  value  of  frequency  ratio  f 

Af  *  increment  used  in  scanning  frequency 

ratios  from  f  to  f. 

a  d 

This  card  contains  the  following  three  dimension¬ 
less  quantities: 

A  =  compressibility  number 

e  “  dimensionless  bearing  steady-state 

x0  displacement 

L/D  “  bearing  length  diameter  ratio 

This  card  contains  the  following  four  control 
integers:  MORE,  MD,  MP,  MACT 


MORE 

-  1, 

another  case  to  follow 

o, 

this  is  the  last  case 

MD 

-  1, 

print  E^,  E2,  E3, 

o, 

print  basic  quantities  only  in 
the  output 

Card  5*  (8E10.3) 

Card  5**  (8E10.3) 


Card  6  (Il.EU.7) 


Card  7 

Cards  6  and  7  may 
Card  6 


MP  =  0,  take  positive  sign  for  complex 

square  root 

-i ,  Lake  negative  sign  for  complex 
square  root 

MACT  -  0,  for  unidirectional  device 

1,  fdr  active  device 

If  MACT  =  1,  two  cards  containing  the  following 
information  are  needed: 


U 

u 

u  , 

u 

mxx 

mxy 

my:: 

tnyy 

V  , 

V 

v 

V 

mxx’ 

mxy 

myx 

myy 

If  MACT  “  0,  one  card  containing  the  following 
four  valueB  is  needed  instead: 

2 

Ag  »  area  of  magnetic  air  gap  [in  ] 

■  area  of  conductor  [cir.  mil] 
l k  ”  length  of  conductor  [in] 
i  ■  current  [amp] 

Change  cards  containing  one  integer  code  and  one 
real  number.  The  real  number  replaces  the  current 
value  of  an  input  quantity  designated  by  the  fol¬ 
lowing  code: 

Code  1 .  i 

8 

Code  2.  C 
Code  3.  a 
Code  4.  e 

xo 

Blank  card 

be  repeated  as  many  times  as  desired. 

An  additional  blank  card  to  terminate 
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The  entire  above  sequence  may  be  repeated  a3  many  times  as  desired  with  MORE  -  l. 


OUTPUT 

Each  output  page  is  for  a  particular  value  of  f,  L'ue  frequency  ratio.  It  will 
begin  with  f  «  f  and  increase  with  increment  Af  until  either  a  critical  frequency 
ratio  is  found  or  the  final  value  f  =  in  reached.  The  following  is  an  illus¬ 
tration  of  a  typical  output  page: 


Line 

1. 

Frequency  ratio 

Line 

2. 

FP  -  F(+) 

FM  -  F(-) 

CAP-F  =  F 

See  Appendix  I 

Line 

3. 

GP  =»  G(+) 

GM  “  G(-) 

See  Appendix  I 

G  -  C 

Line 

4 

K(PARL)  - 

K//" 

Bearing  stiffness  in  line  with  line  of  centers  of 
bearing  and  journal 

F*OMEGA*C ( PASL ) 

“  fa)  C,.  “  Parallel  damping  in  line  with  line  of 
'  centers  of  bearing  and  journal 

K(PERP)  - 

m  Bearing  stiffness  normal  to  line  of  centers  of 
bearing  and  journal 

F*OMEGA*C( PERP ) 

■  f<i)C/  ■  Bearing  damping  normal  to  line  of  cen¬ 
ters  of  bearing  and  journal 

Line 

5. 

W-SUB-XO  - 

W 

xo 

[ dimensionless ] 

W-SUB-YO  =■ 

W 

yo 

[dimensionless] 

ALPHA  *  or 

=  attitude  angle 

F-BAR  *  F 

it 

1 

f  Cw  )2  +  (W  J2]1  ,  [dimensionless] 

1  xo  yo  1  L 

Line 

6. 

F-SUB-M  * 

F 

m 

-  FnTTDLpa,  [lb] 
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K-SUB-M 

U6  1 

II 

-  -2 

C 

F 

m 

[ dimensionless ] 

K-SuB-ma 

kmx 

-  k 

m 

cos 

a 

l dimension less] 

K-SUB-MY 

*  k 

my 

-  k 

m 

•  in 

o 

[ &  1 mens  ion less] 

Lines  7  &  8.  Systems  dynamic  stiffness  and  damping  including  those  contributed 
by  bearing  film  forces  and  magnetic  forces. 

Line  9.  A  * 

8  -  J  B1 

ZR  -  Zr 


Line  10.  W(F)  ■  left  hand  side  of  (11) 

Line  11.  This  additional  line  will  show  only  on^tho  last  output 

page,  if  a  critical  frequency  ratio  is  found. 

M  ■  m  ■  dimensionless  critical  mass 
c 

CAP-M  ■  M  ■  critical  mass  [lb] 

T.S.  ■  threshold  speed  [dimensionless] 

A  Fortran  listing  of  program  PN  424  is  provided  in  the  next  few  pages.  Typical 
listings  of  input  and  output  are  also  given. 
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//  FOR  TOMS 

*I0CS(CARn,14<.3  PRINTER) 

•NUNPROCESS  PROGRAM 
•ONE  WORD  INTEGERS 
♦LIST  ALL 
•PUNCH 
C 

INFLUENCE  OF  MAGNETIC  FORCES 
{UNIDIRECTIONAL  OR  BY  AN  ACTIVE  DEVICE) 
ON  STABILITY  UF  PLAIN  JOURNAL  BEARING 
B.RRYNOLFSON  FOR  T .CHI ANG  2/23/66 


5+  (SE10.3) 


5**  (8610.3) 


INPUT  FOR  T0M5 

1.  <«0H)  TITLE  CARD 

2.  (8E10.3)  SMALL  L *C ,MU, P-SUB-A ,0 

3.  (8610.3)  F-SUB- A , F-SUB-B,  OELTA-F 

4.  (8E10.3)  LAMBDA, EPS  I LON-SUB-XO ,C A P-L/D 

5.  (415)  MURE  *MD»MP, MAC T 

MORE-1  ANOTHER  CASE  TO  FQLiQW  _ _ 

MORE-0  THIS  IS  LAST  CASE 
MO-1  PRINT  E1,E2,E3,E4 
MD-0  ONLY  BASIC  PRINT 

MP-0  +  SIGN  FOR  COMPLEX  SQUARE  ROOT 

MP»— 1  -  SIGN  FOR  COMPLEX  SQUARE  ROOT 

MACT-0  UNIDIRECTIONAL  MAGNETIC  FORCE 
MACT *  1  MAGNETIC  FORCE  FROM  ACTIVE  DEVICE 
IF  MACT* 1 ,  TWO! 2)  CARDS  AS  FOLLOWS 
AUXX «  AUXY ,  AUYX,  AUYY 
A VXX «  AVXY,  A V YX ,  AVYY 
IF  MACT*0,  ONE  CARD  AS  FOLLOWS 
AG,AK, LK, I  _ 

6.  (I1,EI4.7)  CHANGE  CARDS 

7.  BLANK  CARD 

ITEMS  6.  AND  7.  MAY  BE  REPEATED  AS  MANY  TIMES  AS  DESIRED 

8.  AN  ADDITIONAL  BLANK  CARO  TO  TERMINATE 

THE  ENTIRE  ABOVE  SEQUENCE  MAY  BE  REPEATED  AS  MANY  TIMES  AS  OESIRED 
WITH  MORE-1 


CHANGE  CARDS  CONTAIN  ONE  INTEGER  CODE  AND  ONE  REAL  NUMBER. 

THE  REAL  NUMBER  REPLACES  THF  CURRENT  VALUE  OF  AN  INPUT  QUANTITY 
DESIGNATED  BY  THE  INTEGER  CODE. 

THE  CODES  AND  THEIR  C0KRESP0ND1 NG  QUANTITIES  ARE 

1  SMALL  L 

2  C 

3  LAMBDA 

4  EPSILON-SUR-XO 


DIMENSION  W( 30 ) »  FF ( 30 ) 

DATA  NR ,NW/2, 3/ 

DATA  RD2DG,  P I /57 . 29578 , 3. 141 5926535/ 

CALL  TEST 
CALL  DUMPP 
1  READ  (NR, 2000) 


*  .. 


5 


' i7-*K '  'W 


READ  (  NR , 20  10  )  EL,C,ZMU,PSUBA,D 
KFAD  (NR, 2010)  FA,FB,D£LF 

K6AD  (NR, 2010)  ZLAM.tPS.ELOVD  _ 

READ  (NR, 2020)  MORE, HO, MP, MACT 
IF  (MACT)  9,9,5 

5  RFAD  (NR.Prnm  inio.iiiirv.iiivi.iiivv 
READ  (NR, 2010)  A VXX , AVXY , AV YX , AV YY 

go  rn  10 

9  READ  (NR, 2010)  AG, AK, ELK, 01 
10  WRITEINW, 3000) 

WRIT  E  t  NW , 2000 ) 

WRITEINW, 3001) 

WRITEINW, 3010)  EL,C,ZHU, PSUBA,0 

WRITEINW, 3002)  _ 

WRITEINW, 3010)  FA,PB,0EtF 
WRITEINW, 3003) 

WRITEINW, 3010)  ZLAM,EPS, ELOVD 
WRITE (NW, 3004) 

WRITEINW, 3020)  MORE , MO, MP, MACT 
IF  (MACT)  15,15,19 

14  WRITEINW, 3022) 

WRITEINW, 30 10)  AUXX , AUXY , AUYX, AUVY 
WRITEINW, 3023) 

WRITEINW, 3010)  AVXX,AVXY,AVYX,AVYY 
GO  TO  16 

15  WRITEINW, 3025) 

WRITEINW, 3010)  AG, AK, ELK, 01 

16  CARL ■D*ELOVD 
OMEGA-2. 0*C/0 

□MEGA-ZL AM*R SUB A« OMEGA* OMEGA/ ( 6.0-ZHU ) 
WRITEINW, 3017) 

WRITEINW, 3010)  OMEGA 
CALL  ERFUNI ERS ,E1,E2,E3,E4) 

IF  |M0)  22,22,21 

21  WRITEINW, 3005) 

WRITEINW, 3010)  ERS, El, E2.E3.E4 

22  EN-0.0 

I«0  . 

J“1 

25  F«FA*EN*OELF 
FF  I J ) »F 

IF  (FB-F)  150,30,30 

30  WRITEINW, 3101)  F 
CALL  KC  E FG I F  ,  Z  L  AM  ,  E  L  OV 0 ,  F  P ,  FM  , C  A P F  ,GP  ,  GM 

♦FCRRRI 

IF  (MO)  32,32,31 

31  WRITEINW, 30141 
WRITEINW, 3010)  FR,FM,CAPF 
WRITEINW, 3015) 

WRITEINW, 3010)  GP,GM,CAPG 
WRITEINW, 3016) 

WRITEINW, 3010)  ZKRLL,FCPLL,ZKPRP,FCPRP 

32  WXO—0.5*ERS*E1*CARF 
WY0«0,5*EPS*E2*CAPG 
ALFA*ATAN( -WYO/WXO) 

ALF0G-ALFA-RD2DG 
FMBAR-SQRT (WXO*WXO+WYO*WYO) 


WAITE ( N  W , 3006 ) 

WRITE (NW, 3010)  WXO , WYO  *  AL  FDG »  FMBAR 


FM»FMBAR*PI*D*PSUBA*CAPl 

ZKP=-2.0*C*FMBAA/EL 

CSALF«COS(ALFA) 

SNALF-SINI ALFA) 

ZKriX  «Z  Kn»C  S*  i.  r 

zkmy*zkm*snaif  . . . 

CSA2«CSALF*CSALF 

cssn*csalf*$nalf 
SNA2*SNAIF*SNALF 
WRITE (NW, 3007 ) 

WRITE ( NW, 3010 )  FM,ZXM,ZKMX,ZKMY 

UXX  *0 . 5  *  ( E PS*E 3*C A>  fVeF»ZKPL  L)' 

UYX*-0.5*(EPS*E4*CAPG+E2*ZKPRP) 

UXY«0.5*IE1*ZKPRP-EPS*E4*CAPG) 

UYY«0.5*(EPS*E4*CAPF+E2*ZKPLL) 

IF  (HACT)  33,33.34 

MOO  IF  I  ED  4/48  -  NEW ~COEFFTcTENTS  FOR  UNIDIRECTIONAL  CASE 

33  B-S0RTI72.0E4AFM/AG) 

QN*0.62*6*El/QI 
RE*0.B67*0N*EIK/AK 
ELZRO»ON*B*AG*1 .0E-8/0I 

ROVL-RE/ELZRO  __  _ 

WRITE (NW, 3024  f 

WRITE (NW, 3010)  BfONfRE, ELZRG 

FOMEG*F IOMEGA 

t*fomeg*fomeg*rovl*rovl 

FI*  ROVt  *ROVL/T 

F2*F0ME  OVL/T  _ 

WRITE (Nw, 3827) 

WRITE (NW.3010)  F 1 , F2 
T*ZKM*F1 
BUXX«T*CSA2 
BUXY*T*CSSN 

_ BUYX-BUXY  _  _ _ 

BUYY*T *SNA2 
WRITE (NW, 3032! 

WRITE (NW, 3010  I  BUXX,BUXY,BUYX,BUYY 

T*-ZKM*F2 

6VXX*T*CSA2 

BVXY*T*CSSN  _ 

BVYX-BVXY 
BVYY*T*SNA2 
WRITE (NW, 3033! 

WRITE INW.3010 )  BVXX,BVXY,BVYX,BVYY 
UXX*UXX+6UXX 

UXY-UXY+BUXY  _  _  _ 

UYX-UYX+BUYX 
UYY*UYY+BUYY 
WRITE (NW, 3014) 

GO  TO  1035 

34  UXX *UXX+AUXX 
UXY-UXY+AUXY 
UYX  *UYX+AUYX 
UYY*U YY+AUYY 
WRITE (NW, 3021 > 
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1035  WRITE (NW,  3008 ) 

WR ITE (NW, 3010 1  UXX,UXY,UYX,UYY 

C  _ _ _ 

VXX»0.5*El*FCPLl 

VYX=-0.5*E2*FCPRP 

VXY«0.5*E1*FCPRP 

VYY«0.5*E2*FCPLL 

IF  (MACT)  1037,1037,1035 

1035  VXX-VXX+AVXX  _ _ 

VXY«VXY+AVXY 

VYX«VYX*AVYX 

VYY«VYY+AVYY 

GO  TO  1031 
1037  VXX-VXX+BVXX 

vxY«vjdLtevxr  _ 

VYX«VYX*BVYX 
VYY-VYY+BVYY 
1030  WRITE! NW»3009 J 

WRITE (NW, 30 10)  VXX , VX Y, VYX , V YY 
C 

A- ( UXX-UYY  )**l- (  VXX-VYY  ) »*2+4.  0*  ( UX  V»UYX-VXY»V  YX.1 
B«'4.0*(UXY*VYX+UYX*VXY)42.0*(UXX-UYY)*(VXX-VYY) 
CALL  CSORT ( A,B,ZR*ZI ) 

IFIMP+l 134,35,34 
35  ZR--ZR 
ZI--ZI 

_  34  CONTINUE _ _ 

WRITE ( NW, 3010 ) 

WRITE (NW, 3010 )  A,B,ZR,ZI 
W(  JI-VXX+VYY+ZI 
WRITE (NW, 3102 )  W(JI 
C 

JC _ 

IF  (I)  40,40,100 
40  IF  (J-1I  99,99,45 
45  IF  <W(J)*W(J~in  50,100,99 
50  T-IWI  J)-W(  J-1})/<FF(  J)-FFU-1H 
TT«WIJ-1»-T*FF(J-1) 

_ _ 

1-1 

GO  TO  30 
C 

99  EN-EN+l.O 
J«  J+l 

_g. _ _fiQ_I£L25 _ 

C  100  smlm»uxx*vyy+uyy*vxx-uxy*vyx-uyx*vxy 
C  $MLM»SMLH/(F*F*IVXX+VYY ) > 

100  SMLH-  0.5*1 UXX+UYY+ZR I / ( F*FI 
T«C*0MEGA*0MEGA*0 .00259067 

_ EH»PI*CAPL»0*PSUBA»SMLH/T _ 

ts«$qrt<Em*t7?mT 
WRITE (NM,3011 ) 

WRITE (NW, 3010 )  SMLM, EM, TS 
CALL  TEST 
C 

100  CALL  NO R1N(EL,C, ZLAM.EPS, I  ) _ 

IF  (I)  10,10,110 

110  IF  (MORE)  111,111,1 

111  CALL  EXIT 


O’O 


1 50  WRITF!NW,3012>  FA,F8 

_ WR  ITF (NH 1 301 3 )  _ 

J»  J-  1 

00  155  *  1 ,  J 

WHITE (NM, 3030 )  FF(|)tH(l) 
155  CONTINUE 

r n  rn 
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FORMAT S 
2900  FORMAT ( 8 OH 

*  ) 
2010  FORMAT ( IE  10. 3) 

2020  FORMAT (415) 


3000  FORMAT ( /1H1»2*X, ' INFLUENCE  OF  MAGNETIC  FORCES ' ,/21X (UNIOIRECTION 
*AL  OR  BY  AN  ACTIVE  DEV  I CE ) • , /22X , »0N  STABILITY  OF  PLAIN  JOURNAL  BE 
*AR ING* / ) 

3010  FORMAT (1X,8(3XE12«5)) 

3020  FORMAT  (2X, 4(2X15, 2X!5,1X)  ) 


_3.030  .FORMAT!  1* .3X612, 5 ,3ML2.5). 

3001 

3002 

3003 

3004 

3005 

* 


300* 

3007 

I 

3000 

3000 


FORMAT! lHOtlOXIHL, 14X,1HC,1 3X,2HMU»l IX, 7HP-SUB-A.1 IX, 1HD ) 

FORMAT ( 1H0,7X,7HF-SU8-A,8X,7HF-SUB“B,8X,7HDELTA-F> 

FORMAT ( 1H0 *8X ,6HLAMB0A,8X ,7H6PS I  LON, IX  ,  7HCA  P-L/D ) 

FORMAT! 1H0 ,5X (4HM0RE ,  3X , 5HMDI AG* 4X , 2HMP ,4X *4HMACT ) 

FORMAT! 1H0«7X« 

. . _7H.EPS. I  LONtlOX ,2HE1,13X.2HE2,13X,2HE3,13X,2H64) _ 

FORMAT! 1H0.4X, 0HW-SUB-X0,7X,8HW-SWB-Y0, OX, 5HALPHA.10X, 5HF-BAR ) 
FORMAT! 1H0.7X,  > 

7HF-SUB-M.8X, 7HK-SUB-M,7x,8HK-SUB-MX,7X,8HK-SUB-MY) 
FORMAT!  1H0,BX,4HU-XX,11X,4HU-XY,  i;iX,4HU-YX,llX,4HU-YY) 

FORMAT! 1HO,OX,4HV-XX,11X,4HV-XY,11X,4MV-YX, 11X.4HV-YY) 


3011  FORMAT!  1H0,0X  ,1HM.  12X.5HCAP-M,  10X.4HT.S.  ) _ _ 

3012' FORMAT (6H0«**4  ,'W(F)-0  HAS  NO  ROOT  ON  THE  CLOSED  F-INTERVAL  (•« 
•E12.5,1H,,E12.5, 1HJ)  -  * 

3013  FORMAT (lH0f9X,lHF,12X,4HW(F) ) 

3014  FORMAT! 1H0,8X,2HFP,  13X,2MPM,12X,5HCAP-P> 

3015  FORMAT! 1H0, IX, 2HGP,  13X , 2HGM , 14X , 1HG )  , 

301*  FORMAT! 1H0  »*X ,7HK( PARL) «4X,1 5HF»0MEGA»C! PARL)»4X«7HK!PERP).4X, 

*15HF*6mEGA*C!PERP ) » 

3017  FORMAT! 1H0,BX .5H0MEGAI 

3018  FORMAT ( 1H0  ,SX , 1HA , 14X , 1HB  *14X,2HZR*1 3X ,2HZI ) 

3019  FORMAT! 1H0, 'FORCE  COEFFICIENTS  WITH  UNIDIRECTIONAL  MAGNETIC  FORCE* 


*) 

3021  FORMAT! 1H0, • FORCE  COEFFICIENTS  WITH  MAGNETIC  FORCES  GENERATED  BY  A 
*N  ACTIVE  DEVICE  ' ) 

3022  FORMAT! 1H0,*X ,5HAU-XX , 10X , 5HAU-X Y, 1  OX, 5HAU-YX , 10X ,5HAU-YY ) 

3023  FORMAT ( 1H0 ,8X , 5HAV-XX , 10X  ,5HAV-XY, 10X, 5HAV- YX , 10X .5HAV-YY I 
3025  FORMAT ( 1H0, 8X , 2HAG , 1 3X, 2HAK, 1 3X, 2HLK , 14X , 1HI ) 

302*  FORMAT! 1M0 ,9X , 1H6 , 14X , 1HN, 1 3X.2HRE » 1 3X ,2HL0 > 

3027  FORMAT!  1HQ ,BX ,2HF1 , 1 3X, 2HF2 ) _ _ 

3032  FORMAT'!'  1H0  »tX~,5H6U-X"X,l6x  ,  5H6U-XY  ,  10X ,  5HBU-  YX  ,  10X  ,'?HBU-YY  ) 

3033  FORMAT! 1H0.8X ,5HBV-XX , 10X, 5HBV-X Y, 1  OX ,5HBV-YX ,10X , 5HBV-YY > 

C 


3101  FORMAT! /1H1,15X,  3HF  «,E12.5/> 

3102  FORMAT! 1H0 , 12X ,*HW ( F >  -,E12.5/) 

C 

- ENT5 - - — 


i: 


* 

i 

! 

A 


I* 

l 

4 


i 

1 
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//  FOR  MORIN 
♦NONPROCFSS  PROGRAM 

•ONE  WORD  _JNTJLGJELRS  _ _ 

♦LIST  ALL 
♦PUNCH 

SUBROUTINE  MORIN! EL, C.ZLAM.EPS, I ) 

C 

UAIA  NR/2/ 

c  _ _ 

C  CHECK  FIRST  CARD  FOR  BLANK  OR  BAD  INTEGER  COSE 
READ  (NR, 20151  K,VAL 
IF  ! K )  *0,40,1 

1  IF  ( 5-K I  40,40,5 

r. 

C  _  CHECK  SUBSEQUENT  CARO! SI, _ 

2  READ  (NR, 2015)  K, VAL 

IF  ! K )  50,50,3 

3  IF  15-K)  50,50,5  > 

C 

C  IF  O.K.  CODE  GO  TO  MODIFY  CORRESPONDING  INPUT. 

5  GO  TO  (  10, 20, 30. 40).  K _ _ 

lO  EL  -VAL 
GO  TO  2 
20  C  -VAL 
GO  TO  2 
30  ZLAM-VAL 

GO  TO  2  _  _  ~ 

40  EPS  -VAL 
GO  TO  2 
C 

C  END  OF  CHANGE  CARDS,  RETURN  W/  FLAG  SET  TO  RUN  AGAIN. 
50  1*0 
RETURN 

C  FIRST  CARD  BLANK,  END  OF  CASE! 

40  1*1 
RETURN 
C 

2015  FORMAT ( 1 1 ,E 14. 7 ) 

ENO  _ _ _ 


//  FOR  EPFUN 
♦NONPROCESS  PROGRAM 


♦ONE  WORD  INTEGERS _ 

♦LIST  ALL 
♦PUNCH 

SUBROUTINE  EPFUN ( E , El , E2, E 3 , E* ) 
IF ( E )  1,1,2 

1  fcl*l.tJ 

_ L2ii^ _ 

E3-0.0 
EA-fl.O 
GO  TO  3 

2  EE*E*E 
EEE-EE*E 

EH«1.0-ES _ 

ES-SORT  (EM) 
tS3«ES*EM 
E2«2.0/EE*( l.O-ES) 

E1-E2/ES 

EA*( 1 »0-ES3)*2.0 

£S«2»Q/EEE/ES3 _ 

E3-E5*(3.0<*EE-E«) 

E4«E5*  <  EE* ( EE -3 .0 ) +E4 ) 

3  RETURN 
ENO 
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//  fur  Kcer-G 
•NONPRUCESS  PROGRAM 

•  ONE  WORD  INTEGERS  _ _ _ 

•LIST  ALL 
•PUNCH 

SUBROUTINE  KCEFG) ER ♦FLAM, ALOD,FP, EM, PO ,GP ,GM , GO, AK1  , AC  I , AK2, AC2 ) 

w  _  * 

t\—  -  A 

PcPLAM 

3  P2-1.0-P-P  _ _ _ _ 

PS-S0RT  (P2) 

PF«  P/P2 
AL*0  #5* ( PS+1 ,G ) 

BE*0.5*(P$-1.0) 

BE-  SORT  (BE  I 

AL-SOBT  (ALj _ 

IF ( P  I  1,1,2 

1  BE— BE 

2  ABP*AL-B£»P 
APR=AL*P+SE 

IF ( ALOD- 100.0  I  T,S,B 

___-7  PS*ALOD»PS  _  ______  _ _ 

ARG1*2.0*AL*AL0D 
ARG2«2.0*BE*AL0D 
E-EXP  ( ARG1 ) 

SH-0.5*(E-1.0/E) 

CH*0,5*( E+l.O/E  ) 

SF-SIN  (ARG2I  _  _ 

CF-COS  <  ARG 2 ) 

DEN-P  S*  <  CH+CF ) 

a  XF-PF-P 

XG— PF 

IF ( AL 00-100. 0 )  *,10,10 

?  XF«PF*(ABP*SF-APB*SM)/DEN±_XF _ 

XG-  PF*< ABP*SH+APB-SF)/DEN-XG 
10  1F(K)  4,5,* 

4  FO-XF 
GO*-XG 
.K-0 

_  ej»Pl***<l*0-42.0*FRJ[ _ 

GO  TO  3 

5  FP*XF 
GP*-XG 
K«1 

P-PL  AM* ( 1 ,0-2»0*FR I 

GO  TO  3  __  _ 

A  FH-XF 
GM--XG 

AK1«0.5*(FM+FP) 

AC  1  —0 •  5  * ( GM-GP ) 

AX2»0.5*(GM+GP) 

AC2«0.5*(FM-FP) 

. . RETURN  ■  '  . . 

END 


'  3 
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i :  * 
.  i. 


INFLUENCE  Of  MAGNETIC  FOPCES 
UIN!  PT  PE  CT  TONAL  PP  ?T  AN  ACTIVE  OE  VI  CE_I_ 
nr  STABILITY  0*  PLAIN  JOUPNAl  BEAPING 


sample  p&  pp  i  ir  t  t  on 

PUN  JUNF 

l  <jf  g 

i 

C 

MU 

P-SU B- A 

0 

.  i  cp  on  -r>  i 
r  -s  yp  -a 

. 20  CO  0-02 
F-SUB-B 

.2  70  00-0  8 
1  tL  TA  -f 

.IN  70  04  02 

•  3  SO  00  +0  1 

.  3  6p  m  ~n  n 

.  «*c  CO  C-  00 

.2  00  00  -0  l 

L  A  mb  o* 

CPSILON 

C  A»  -L  /D 

.1  n?  nD  *o  i 

.  3«  m  n-on 

.1  00  00  *0  1 

KOBE  HP  lie 

1  o 

HP  MAC  T 

-l  l 

tri  -*  x 
.2  TIP  CP 

Au-  XV 

.no  ooo 

»u-vx 
.0  00  00 

A  U-  VV 

.00  00  0 

~  nr-rr 
.  *  oc  on  -n  i 

rv-ntv 
.no  ooo 

AV-TV - 

.0  00  00 

- imrT 

.  30  00  0-01 

CfTP* 

.1  2S  Aj  *0  1, 

- 

t  . 

t 


I  i 
■ 

i  ' 


i  y°  jff  *—'*-• 


r 


3«  am- oo 


W-  Slip-  xo 

w-sup-yo 

AL  °H  A 

F-BAT 

-.1  54  33  -0  1 

.46  38  1-01 

.7 18  IS ♦0  2 

.  49  45  1-  01 

r  -«•  on  -»* 

K- V,MI-M 

N-  SI  IP  -  MX 

K  -XllB-M  y 

.2  79  76  .p? 

-.19  90-ni 

- .6  17  33-n  7 

ift  7Q  m 

F0.7CF  COFFFTCIFNTC 

V  IT H  MAGNET  IF 

FflOCES  CENfpATEO 

BX  AN  ACTIVE  DEVICE- 

li-XX 

u-x  x 

u-  rx 

t'-V  y 

,6  39  irt-ai 

•  10  52  8*00 

-  .1  16  79*00 

.52  45  3-  01 

X-  XX 

tf-X  V 

y-  xx 

V-Y  y 

•  i  ii  i«  *m 

-.si  (71  9-01 

.4  70  39-0  1 

.  in  48  5*  00 

A 

Q 

2» 

2  T 

-  .  3  9*1  <13  -o  ! 

.93  786-ai 

-.986  72-0  1 

-.2218  8-00 

WlF)  t 

-.50 96 5- 02 

F  z 

.4  0  00  o-on 

v-sun-  xn 
-.1  54  33  -Hl 

W-SllR-YO 
.46  98  1-01 

ALPHA 

.71815*32 

r-a  ap 

.49  45  1-01 

F  -S  08  — M 
.2  79  76*0  2 

K-S1J9-M 
-.19  78  0-  01 

K-SUh-MX 
-.6  1733-0  2 

K-StIB  -MY 

18  79  2-  01 

FOPCF  CPEFP!  'TFNTS 

WITH  MAGNETIC 

FOPCES  GENFPATEO 

GY  AN  ACTIVE  OEVICE 

U-  XX 

.6  5369-01 

crnry 

.10  36  6*  00 

IKYX 

-.1  IS  30*00 

U-YT 

.53  79  7-01 

V-  XX 

.1  15  24  *00 

V-X  V 

-.53233-01 

V-YX 

.4  90  85  -0  l 

V-Y  Y 

.  10  85  9*  00 

4 

-.3  72  64  -01 

B 

.45  056-01 

29 

-.10297  *00 

ZT 

-.  21  87  8-  00 

BCF1  r 

.50456-02 

APPENDIX  IV:  STABILIZATION  WITH  QUADRATIC  DAMPING 
SPECIAL  SYMBOLS  FOR  APPENDIX  IV 

Noro:  The  unit  is  dimensionless  unless  otherwise  indicated  in  brackets. 


A; 

& 


f 

L 

M 


A 


xxn 


radial  bearing  clearance  [in.) 
bearing  diameter  [in.) 

dimensionless  work  per  cycle  dissipated  by  virtual  damping 

Cl/ Ui,  frequency  ratio 

bearing  length  [in.1 

rotor  mass  [ lb- sec  /in. “1 
2 

HCo)  /(p^nLD) ,  dimensionless  rotor  mass 
harmonic  number  of  the  shaft  center  orbit 
ambient  pressure  [pais) 
time  [sec.) 

peak  amplitude  of  the  in-phase  component  of  W  in  response  to  the  nth 


xyn 

I  u 

)  yxn 

luyy«. 

U.n 

^n  J 


,  u 

/  xn 


v 

xn  J 


harmonic  of  the 


U 


(!) 

»f  the 

of  ch'(‘) 


component  of  the  shaft  center  orbit 


peak  amplitude  of  the  in-phase  component  of  W  .  in  response  to  the  nth 
f  > I  yt 

harmonic 

r 


component  of  the  shaft  center  orbit 


xxn 


xxn 


yyn 


yyn 


U  -  -u 

xyn  yxn 


for  e 


xo 


yn 


yn 


xyn  yxn  J 

Fourier  coefficient  of 

Fourier  coefficient  of 


cos  nft  | 

sin  nfT J  ^Wxt^m 

cos  nfT""^ 

J  >  of  (W  J 

sin  nf t  f  yt'm 
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to 


yn 


n 


phase  angle  of  the  nth  harmonic  of  the  y  component  of  the  shaft 
center  orbit  [radians] 

fundamental  orbit  frequency  of  shaft  center  [radians/sec .  J 
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The  stability  analysis  cited  in  Section  II  and  Appendix  I  is  applicable  only  to 
linear  systems,  and  therefore  it  is  not  capable  of  coping  with  the  concept  pro¬ 
posed  in  Section  III  B,  in  which  electro-magnets  would  be  used  to  impcjse  quadratic 
damping  forces  on  the  shaft.  In  this  Appendix,  a  more  general  criterion  to  (^- 
termine  the  stability  of  the  shaft-bearing  system  will  be  derived.  The  more.? 
general  criterion  is  based  on  the  concept  of  "virtual  damping  for  a  stationary 
state" 

Consider  the  dynamic  equilibrium  amongst  the  virtual  damping  force,  the  D'alembert 
force  of  the  shaft,  and  the  bearing  force  end  the  force  of  the  electro-magnet  in 
the  dimensionless  form: 


T  *  ”«  +  ■ 


(IV.  U) 


*  “yt  ♦  <V« 


(IV.  lb) 


0  is  the  "virtual  damping  coefficient"  (written  non-dimensionally) .  It  Is  assumed, 
that  with  the  introduction  of  the  virtual  damping,  the  system  would  sustain  a  state 
of  periodic  motion  such  that 


«x  (T  *  2rr/f)  =  (t  -  0) 


€y  (e  *  2ir/f)  -  *y  (T  -  0) 


(IV. 2) 


The  energy  dissipated  by  the  virtual  damping  per  cycle  of  the  periodic  motion  is 
obtained  by  multiplying  d*  /dT  and  da  /dT  into  Eqs.  (IV. la)  and  (IV. lb)  then  in- 

*  y 

tegrating  over  one  period: 


2rr/f 

$•  \  • 

0 

f  de 

1  ,• — 

2 

+ 

de  2~| 

J 

dT 

H 

• 

a 

00-0 

N. 

d2eK 

(dr2 

de 

X 

d2e 

de 
_ 1 

dT 

dT2 

dT  ‘ 

+  'i  k"xt+  ‘V.fe5'  +  (v  dT  I 

o  L 

d  e  de  .  dc  *  d*e  de  ,  de  £ 

Since  — -"  dT  *  —  d  (t— )  and  — 77^  dr  ■  -  d  (t“)  and  since  the 

.2  dT  2  dT  ,2  dT  i  dT 


d«x 

dc 

periodic  condition  applies  to 

and  -r-^  as  well 
aT 

d2,  d«  d2. 

dc 

•  s  # * + 

d^>  dT  "  0 

Clearly,  this  simply  indicates  that  the  D'alembert  forces  would  alternately 
store  and  release  energy  during  different  parts  of  the  period,  but  would  have 
no  net  contribution  to  the  energy  of  the  system  when  the  entire  period  is  con¬ 
sidered.  Thus 
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If  A  is  zero,  the  prescribed  motion  (e  ,  s  )  defines  a  stationary  state  of 

x  y 

periodic  motion.  If  0  is  positive,  implying  the  need  of  additional  damping  to 
sustain  the  prescribed  motion,  the  system  Is  unstable;  and  if  0  is  negative, 
conversely,  the  system  is  stable. 

In  general,  one  must  permit  the  periodic  motion  to  contain  harmonic  components, 
e.gs 


«x  “  ‘xo  +  2  *xtn  Cofl  (n£T  ’  V 


I  <  „  cos  (nfr  -  a  ) 
ytn  yn 


(IV.  4) 


Coupling  between  the  two  degrees  of  freedom,  x  and  y,  is  usually  Imposed  through 
the  bearing  forces,  Wxt  and  Wyt>  both  of  which  generally  depend  on  both  and  *^. 


Assuming  that  c  and  c  are  aaall  enough  so  that  the  orbit  amplitude  is 

...  1 1  &  Jr  ••  1 1 

small  in  comparison  with  the  bearing  film  thickness,  then,  consistent  with 
Equation  (IV. 4) ,  the  dynamic  bearing  forces  can  be  expressed  as: 


1 

xxn 

cos 

(nfT  - 

®  ) 

Txn 

-  V 

xxn 

1 

xyn 

cos 

(nfT  - 

'P  ) 

- 

yn' 

xyn 

1 

yxn 

cos 

(nfT  - 

cp  ) 

-  V 

xn 

yxn 

I 

yyn 

cos 

(nfT  - 

<p  > 

-  V 

yn 

yyn 

J  eytn]j 


(IV. 5) 
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are  solutions  of  the 


For  an  example,  consider 


e  “0 
xo 


(IV. 8) 


consequently,  due  to  symmetry, 


u 

xxn 

-  U 

yyn 

m 

n 

V 

xxn 

=  V 

yyn 

m 

% 

U 

xyn 

*  -U 

yxn 

SB 

Uin 

V 

xyn 

-  -V 

yxn 

at 

V^n 

(IV. 9) 


which  are  dynamic  perturbation  aolutiona  corresponding  to  the  frequency  of  the 
unloaded  plain  journal  bearing  cited  in  Appendix  I  and  [Ref.  5].  Thus, 


Wxt  *  *  0//1  cos  fT  ‘  V//l  8ln  fT)]‘xtl 

*  [<^1  cOS  <fT  *  ®yi>  ■  ^1  sln  (fT  "  ®yl0  ®ytl 

Wyt  *  C  Ual  cos  fT  *  VUl  sln  fT3*xtl 

-  [o#l  cos  (fT  -  oyl)  -  V/n  sin  (fT  -  ®yl)]  eytl 


(IV. 10) 


Differentiating  Equations  (IV. 7)  with  respect  to  t: 


de 

*'  -  f  «  sin  fT 
dT  xtl 

de 

d^  '  -  f  *ytl  8in  (fT  •  V 


(IV. 11) 
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.*4  **■•#■*' 


To  make  3  most  positive,  one  finds 


£Xtl  "  6ytl  “  61 


(IV. 16) 


and 


sin®yj  -  Bg(uul) 


Actually,  since  is  always  positive  according  to  numerical  results,  one  simply 

has 


V  2 


:  <jtV.l7) 


Thus . 


V//l  +  Uil]  ‘  3rr  em  f  *1 


{IV.  18) 


At  the  stationary  state,  8*0,  then 


(«,) 


1  stationary  6 


(H/1  ~  v//l* 


KS 

m 


(IV. 19) 


If  <  ^ei)stationary>  P  would  be  positive,  then  the  system  is  unstable,  and 

would  grow.  If  >  (,i)sta£ioIiary  P  would  be  negative,  and  «j  would 

diminish-  Thus  the  presence  of  (W  .)  and  (W  .)  stabilizes  the  orbit  sice  at 

xt  m  yt  m 

€j,  which,  however,  can  not  be  completely  reduced  to  eero.  To  determine  f,  sub¬ 
stitute  Equations  (IV. 6),  (IV. 7),  (IV. TO),  (IV. 11),  (IV. 16),  (IV.17)  and  8  -  0  into 


Equation  {.IV.lalr 


'-'/I 


+  V^  sin  ft  -  Uj ^  sin  ft  -  V, j  cos  £t 


+  P  f  e,  sin  £t 

iP  1 


sin  fT 


si 


bi,[' 


cos  fr 


“u  -  V  ■ e. £  S1 


“  -\ 

%  *  VJ 1  -  "f2 

—  w 


sin  ft 


j 


(IV. 20) 


One  should  note  that  the  last  expression  cannot  be  precisely  zero  because  the 
truncated  description  of  the  periodic  motion  is  only  an  approximation.  Also, 
Eq.  (IV  .lb)  needs  not  be  separately  considered  becaulse  of  the  prevailing  sym¬ 
metry.  Since,  even  and  odd  functions  of  t  should  separately  vanish  in  the 
above  expression,  one  would  find 


U//l  +  V^1 
m 


(IV  21.) 


Harmonic  contents  in  the  motion  of  the  stationary  state  can  be  obtained  from 
Eqs  (IV.  3.  a,b)  by  performing  a  Fourier  analysis  in  an  iterative  manner  begin¬ 
ning  with  the  truncated  expressions,  Eq.  (IV. 7),  as  the  initial  guess  in  the 
non-linear  terms,  which  are  defined  by  Eq.  (IV. 6).  The  iterative  procedure  may 
be  continued  to  improve  accuracy  further-  The  process  begins  with  finding  the 
harmonic  concents  of  Eq .  (IV. 6)  by  the  substitution  of  Eqs.  (IV. 11),  (IV. 16) 
and  (IV. 17)  and  then  performing  the  appropriate  Fourier  analysis; 
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(u 


xt  m 


n=  1 


xn  cos  nfT  +  sin  nfT) 


(w...)_a  (u 


yt'm 


n«*l 


yn  COS  nfT  +  vyn  Sin  nfT) 


(IV. 22) 


2n/f 


a  i  c 

m  T7  J 

0 


l!sl 

dT  1  dT~  C°®  nfT  dT 


2rr/f 


pi  C  |!i 

”  "  J  d  t 

0 


xn 


dc 


37“  sin  nfT  dT 


yn 


A 

3  i  „C 

m  tt  J 


2n/f 


!  d«  d« 

V  I  _  y 

i  dT  '  dT  C0B  nft  dT 


yn 


3  -  [ 

®  TT  J 


2n/f 


d* 

—2 

dT 


dc 

37*  sin  nf T  dT 


1 


)  (IV. 23) 


J 


0r}  with  the  aid  of  Eqs.  (IV.1I)  and  (IV.17) 
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ri 


Substituting  Eqs.  (IV  4),  (IV. 5),  (Iv-9).  and  (T.V.22)  intoEq.  UV.I/ii  (iy.lt) 
netting  p  -  l)  for  the  stationary  state - 


r-  2  2 

m  ,  n  f  e  .  cos  (nfT  -  ©  ) 

i-  xtn  xn 

n-1 


at 

^  f  e  fu 
L.  |  xtn  ^  //n 

n»i 


cos 


<nfT  -  raxn)  -  V//a  sin  (nfr  -  ®xn)] 


xn 


+  *ytn  COS  (nfT  *  <0'  '  V'~  eln  (nfT 


yn  'in 


•  V^} 


v 


+  ,  (u„  c09  nfT  +  v  sin  nfT) 

i~  xn  xn 

n*l 


V 

n»l 


(cos  nfT  fOnf  .sin  ©  -  cu  n2f2  cos  +  U„  cos  ©  +  V.,  sin  eg  )t 

u  xn  xn  //n  xn  //n  xn' 


+  (U,  cos  ©  +  V.  sin  <g  )  c  „  -  u  ™| 

-fn  yn  -*n  yn  ytn  xnj 


+  sin  nfT  [(-Bnf  cos  ®xn  -  mn2f2  sin  ©xn  +  V  sin  <0  -  V//n  cos  ©  )c, 


xn  xtn 


+  (Vj  sin  co  -  V.  cos  ©  )  * 

-*n  yn  4  n  yn  ytn 


v  “A 

xn  J  I 


xtn 
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Or, 


IT  1 

//n- 


ms  m 


U 


2,2 

m  n  f 


"xn 


if  **  *•  * 

7/n  '“xn  J ’xtn  - 


-  (U  cos  co  +  V,  sin  co  )  €  „ 
in  yn  in  yn  ytn 


xn 


>  (IV. 25a) 


f(m  n2f2  -  U„  )  sin  to  +  V,,  cos  co  He 
*-  //n  xn  //n  xn  J  : 


xtn 


(U  sin  co  -  V  cos  co  )  c  _ 
jn  yn  in  yn  ytn 


xn 


and  similarly. 


,  2.2 

(m  n  f 


W  C0B  > 


V//n 


sin  co 

yn 


<D2n  C0S 


CD  + 

xn 


V^n  8in  ®xn> 


xtn 


i-  u 


yn 


)  (IV. 25b) 


P (m  n2f2  -  l  )  sin  co  +  V..  cos  co  "1  * 
//n  yn  //n  ynj 


ynj  ytn 


+  (U  sin  co  -  V,  cos  co  )  c  . 
in  xn  in  xn  xtn 


>-  v 


yn 


It  is  not  necessary  to  consider  Eqs.  (IV. 25)  with  n  -  1  since  these  conditions 
would  be  consistent  with  Eqs.  (IV. 7),  (IV. 16),  (IV. 17),  and  (IV. 19).  For  n  =  an 
even  integer, because  uxn“. Uyn~  vxn  *  vyn  ■  0,  one  also  must  have  extn  «*  eytn~  0. 

For  n  -  an  odd  integer,  Eqs.  (IV. 25)  can  be  rewritten  into  the  following  matrix 
form: 
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'1 


*  tan 


■i  (°  -2f2  •  »//n  -  y 

•  V 


(IV. 27) 


and  by  virtue  of  Eq.  (IV. 19), 


e  *  c 
xtn  ytn 


n(4-n^) 


<".11  •  W 

(m  U  -  V  )*1  +  (U  V  7 

in  .in  //n 


(IV. 28) 


These  two  formulae  are  valid  for  all  positive  odd  integers,  including  n  *  1.  In 
principle,  additional  Iterations  can  be  performed. 

Summarizing,  according  to  the  one-step  iteration  analysis  carried  out  above,  an 
unloaded,  plain  journal  bearing-shaft  system  would  be  stabilized  by  the  electro¬ 
magnets  to  assume  a  "steady-state"  orbit.  Collecting  the  relevant  formulae,  Eqs. 
(IV. 21),  (IV. 19),  (IV. 27),  and  (IV. 28)  and  rearranging  somewhat  for  clearer  pres¬ 
entation,  the  "steady-state"  orbit  can  be  described  as  follows: 


3tt  "fl  '  fyl.  .■  r 

e‘v  +  v*„  A 


n-1,3,5' 


n(4-n^) 


<”il  -  V^)  cos  (nfT  -  con) 

f(m  n2f2  -  U  -  V  )Z  +  (U  -  V. 

//n  xn  j.n  //a 


3rr  Uil  '  V//l.  _m  T 

8  U//l  +  -V  8m  L 


L 

n-1,3,5' 


n(4-n^) 


»"f2  •  -  vi„>J  +  <»,„  -  vA 
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(IV. 29) 
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